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Résumé 


L'objet général de notre travail est l’étude de l’existence globale -en temps- des 
solutions pour des systèmes de Réaction-Diffusion à réaction fractionnaire de 
type Gierer-Meinhardt. 

Ce type de systèmes a une vaste application aux domaines de la Biologie et de la 
Biochimie. 

Nous avons étudié l'existence globale des solutions pour un certain système R- 
D de type Gierer-Meinhardt composé de trois équations. 

Par ailleurs, nous avons généralisé l'étude de l'explosion des solutions pour le 
système EDO associé au système de Gierer-Meinhardt à m équations. Ce modèle 
est composé d'un activateur et de (m-1) inhibiteurs. 

Nos techniques sont basées sur la construction de fonctionnelles de Lyapunov, 
les semi-groupes, les formes quadratiques, etc. 

Mots-clés: Systèmes de Réaction-Diffusion, systèmes de Gierer-Meinhardt, 
modèle activateur-inhibiteur, existence globale, explosion en temps fini, 
fonctionnelle de Lyapunov, semi- groupe. 
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Abstract 


The general purpose of our work is to study the global existence -in time- of 
solutions for Reaction-Diffusion Systems with fractional reaction of Gierer- 
Meinhardt type. 

This type of Systems has a wide applications in the fields of Biology and 
Biochemistry. 

We studied the global existence of solutions for a System of Gierer-Meinhardt 
type of three équations. 

Next, the study of blow up of solutions was generalized for the ODE associated 
with the Gierer-Meinhardt System correspond of m équations. This model is 
composed with one activator and ( m-1 ) inhibitors. 

Our techniques are based on the construction of Lyapunov functional, semi- 
groups, quadratic forms, etc. 

Keywords: Reaction-Diffusion Systems, Gierer-Meinhardt Systems, Activator- 
Inhibitor model, Lyapunov functional, semigroup, global existence, blow up. 
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0.1 Introduction Générale 


Les systèmes de réaction-diffusion sont des systèmes d’équations différentielles couplés 
aux dérivées partielles de type parabolique. Aujourd’hui, ces systèmes reçoivent un 
grand traité d’attention, pour la richesse de la structure de leurs ensembles de solu¬ 
tions. De tels systèmes modélisent des phénomènes apparaissant dans des secteurs 
variés: chimie, biologie, neurophysiologie, épidémiologie, combustion, génétique des 
populations, etc. Les individus diffèrent d’un problème à un autre: en chimie, par ex¬ 
emple, ils représentent des substances chimiques; en biochimie, ils peuvent représenter 
des molécules; en métallurgie, des atomes; en dynamique des populations, ce sont des 
humains; en génétique des populations, ils représentent des caractères. En biophysique, 
ils représentent des charges électriques ou bien des différences de potentiel et en envi¬ 
ronnement, ce sont des animaux ou des plantes d’une forêt, d’une mer ou bien d’un 
fleuve. 

L’étude du présent travail est l’existence globale et/ou l’explosion (blow up) en 
temps fini des solutions de Systèmes de Réaction-Diffusion de type Gierer-Meinhardt 
formés d’équations aux dérivées partielles de type parabolique avec une réaction frac¬ 
tionnaire. On utilise dans cette optique quelques méthodes fonctionnelles et algébriques 
telles que la construction des fonctionnelles de Lyapunov, les semi-groupes, les formes 
quadratiques, etc. 

La thèse est structurée comme suit: 

Chapitre. 1: On rappelle les notations et les notions préliminaires nécessaires 
en vue de les utiliser dans les chapitres suivants. On insistera, en particulier, sur 
quelques théorèmes généraux d’analyse fonctionnelle, sur les espaces de Sobolev, les 
semi-groupes, quelques inégalités utiles et la positivité des formes quadratiques. 

Chapitre 2: On donne des applications concrètes des systèmes de réaction-diffusion 
dans plusieurs domaines du monde réel. 

Ainsi exposera des méthodes générales transformant les phénomènes naturelles à des 
Systèmes de Réaction-Diffusion de la forme 

(|-DA« = /(«)), 
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en utilisant la seconde loi de Fick. 

Par ailleurs: On étudiera quelques exemples dans des domaines variés: En bioécologie, 
électronique et en génétique des populations. 

Enfin, on donnera un exemple d’étude de l’existence globale en temps des solutions 
d’un système de réaction-diffusion avec une matrice de diffusion triangulaire. 

Chapitre 3: Dans ce chapitre nous traiterons un système de Gierer-Meinhardt qui 
est un type célèbre parmi les systèmes de Réaction-Diffusion. On donnera ensuite une 
forme pratique de ce système modélisant un certain phénomène biologique. D’autre 
part, on présentera un exemple de modélisation pour la coloration des coquilles de 
mer. Finalement, on étudiera l’existence globale des solutions d’un système de Gierer- 
Meinhardt à deux équations. 

Au Chapitre 4: Nous établirons l’existence globale des solutions d’un certain 
système de Réaction-Diffusion de type Gierer-Meinhardt composé de trois équations, 
appliqué à la biologie végétale; au développement des plantes dans un domaine appelé 
phyllotaxie (disposition des palmes sur le stipe de la plante). Pour cette raison on 
appliquera quelques lemmes auxiliaires d’estimation ainsi que des théorèmes qu’on 
prouvera à l’aide de la construction d’une fonctionnelle de Lyapunov. Nous allons 
également discuter l’explosion -en temps fini- du système EDO correspondant. 

D’autre part on étudiera l’existence globale des solutions d’un système couplé à quatre 
équations de type Gierer-Meinhardt modifié à l’aide de la modification de Conway- 
Cooper sur la troisième équation. 

Finalement, dans le Cinquième Chapitre, on généralisera l’étude de l’explosion 
des solutions du système EDO associé au système de type Gierer-Meinhardt à m équa¬ 
tions ou encore d’un modèle composé d’un activateur et de (m — 1) inhibiteurs. 
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Chapitre 1 


Notations et Notions Générales 


1.1 Notations générales 


0 étant un domaine borné non vide de WL N , de frontière <90. 
da désignera la mesure superficielle sur <90. 

On définit le gradient d’une fonction u .... x n ) par: 


gradw = Vu 


du du du 


et son Laplacien par 


Au = Y.] 

i= 1 


d 2 u 

dxf 


( 1 . 1 . 1 ) 


— est la dérivée normale de u extérieure à <90. c’est-à-dire: 
drj 

du _> 

^ = v “' "■ 

où Tf désigne le vecteur unitaire de la normale extérieure à <90. 

(7(0) est l’espace vectoriel des fonctions continues dans O. 

C fc (0) est l’espace des fonctions k fois continûment différentiables sur O, où k est un 
entier positif non nul et on a 


C°°(0) = p|C fc (0). 

k> o 
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Un espace de Banach est un espace vectoriel normé complet i.e. toute suite de 
Cauchy est convergente. 

S)(O) est l’espace des fonctions C°° à support compact contenu dans Q. 

S)'(O) est l’espace des distributions sur Q c’est-à-dire l’espace des formes linéaires 
continues sur 5D (C2). 


1.1.1 Espaces fonctionnels 


L p (f2) est l’espace des fonctions u mesurables sur 0 telles que f \u\ p dx < oo, 1 < p < oo 

n 

muni de la norme 


u 


i p - 


I LP (fl) 


= W\Ja dX ■ 


( 1 . 1 . 2 ) 


L°°(0) est l’espace des fonctions u mesurables telles que |tt| < M, p.p (presque partout) 
sur O, où M est une constante. L°°(0) est muni de la norme 


|'u||loo ( -q- ) = supess |u (x)| 


(1.1.3) 


On a aussi, si (X, ||.|| x ) est un Banach; 

L p (0,T, X) = | ri mesurable de [0, T]—> X, fj' || u\\^dt < oo,si 1 < p < oo j muni de 
la norme: 


u 


LP(0,T,X) 


IwIIy dt 


(1.1.4) 


L°°(0,T,X) — lu mesurable de [0,T] —> X, supess ||«|| x < oo > muni de la norme: 

te(o,T) 


\ u \\l~(o,t,x) = supess II u\ 
te(o,T) 


(1.1.5) 


C 0,a (U) = <| u G C (U) ; sup < oo \ avec 0 < a < 1. 


x,y€Q 


c k ' a (U) = {uec k (U) ; G c°’ a (U) Vj, \j\ < k) . 

On définit l’espace H 1 (U) par 

r du 

tf(0)= MG L 2 (U), — G L 2 (U), 1 < i < n 

I •Xj q 
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muni de la norme 


u 


Ih 1 ^) 


= / \u\ dx + 


r» n 

E 

i— 1 


du 


d Xi 


dx = 


u 


' dx + / | Vu| 2 dx. 


( 1 . 1 . 6 ) 


Généralement on a pour m G N* 


H™'(O) = [u G IL 2 (O), D a u G IL 2 (O) pour tout a G vérifiant |a| < m} 
muni de la norme 

IMlH m (Q) = 

\a\<m Q 

= E ll B ““fe(«) • (1-1-7) 

|a|<m 

^ai+û!2 + '”+û: n n 

où ü“ = - „ „ —--et lai = Y )ati la dérivée au sens des distributions. 

dx^dx% ...ox% n ~i 

Définissons l’éspace W m - p (0) (1 < p < oo) par 


W m ' p (0) = {u G 1/(0), D a u G 1/(0), |a| < m} 
muni de la norme 

NI m.p= ll^llu-cn) 1 l<P<+oo- (1-1-8) 

\a\<m 

On rappellera que H 1 (O), H" 1 (O), W m,p (O) sont appelés des espaces de Sobolev. 


1.2 Notions générales 

Déterminants principaux successifs d’une matrice 
Soit la matrice carrée d’ordre n 


( 


( a i,j)l<i,j<n A 


®11 

®12 

Clin 

0-21 

O 22 

a 2 n 

Onl 

On2 

&nn 


\ 


(1.2.9) 
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Ses déterminants principaux successifs sont notés par detl, det2,..., detn, et 
écrits sous la forme: 


det 1 = an, det 2 = 


an «12 
«21 «22 


,det n = 


a u 

«12 

Clin 

«21 

«22 

n ,, _ 



(1.2.10) 

«ni 

«n 2 



1.2.1 Formules de Green 

Nous rappelons maintenant quelques formules de Green qui généralisent au cas multi¬ 
dimensionnel la formule d’intégration par parties en dimension une. Elles sénoncent- 
de la manière suivante: 

Théorème 1.2.1 On suppose que fl est un domaine ouvert de frontière d fl continue par 
morceaux. 

Alors, si u et v sont des fonctions de H 1 (fl) , on a 


'du f dv f 

——vdx = — u——dx + / uvn-da , 1 <i <n, 
dxi J dxi J 

n an 


( 1 . 2 . 11 ) 


fl dfl 

on désigne par g, le i eme cosinus directeur de la normale g à dfl dirigée vers l’extérieur 
de fl et on écrit g i = ( rpel ), da la mesure superficielle sur dfl. 

Preuve. Si u (resp. v ) appartient à H 1 (fl), il existe une suite ( u m ) (resp. (v p )) 
de S)(fl) qui converge vers u dans H 1 (fl) (resp. vers v dans H 1 (fl)) [D(fl) dense dans 
IHI 1 (fl)]. 

On a pour les fonctions u rn et v p de 2) (fl) 


J ~^x^ Vp< ^ X = ~I U ^ X + J 1 - * - n -> 

n n an 


( 1 . 2 . 12 ) 


on obtient l’expression (1.2.11) par passage à la limite dans la formule précédente. 

Corollaire 1.2.1 Pour toute fonction u de EI 2 (fl) et toute fonction v de EI 1 (fl), on a la 
formule de Green 

J (A u) v = j —vdcr — J VaVn. (1.2.13) 

fl dfl fl 



Preuve. Donnons une conséquence de Théorème (1.2.1). 

n d 2 u 

On pose Au = y'—— le Laplacien d’une distribution u. Alors, si u est une fonction 

i=i dx i 

de H 2 (f2), on a d’après (1.2.11) pour toute fonction v de H 1 (f2) 


-j(Au)v = -±j 
n %=r ' 

n 

= E 

i —1 
n 

= E 


d 2 u 

dx 2 


vdx 


du du 


dxi dxi 


dx — [ ——rijjda 


dxi 


n 

du du 


an 


i =1 ' 


dxi dxi 


dx- f^vda 


an 


drj 


'du 


= / VuWvdx — I —vdcr. 

drj 


an 


Remarque 1.2.1 La formule (1.2.11) reste valable si u,v G C l et la formule (1.2.13) 
reste valable si u E C 2 ,v G C 1 . 

Remarque 1.2.2 les formules de Green sont des généralisations de la formule d’intégration 
par parties de la dimension une à la dimension n. 


1.2.2 Seconde loi de Fick 

Considérons le flux de particules d’une certaine espèce, les particules peuvent être 
des molécules, des atomes, des défauts ponctuels, des électrons libres ou des trous 
électroniques, etc. Soit C(x, t ) leur concentration, exprimée en nombre de particules 
ou d’atomes par unité de volume. 

On définit le flux de diffusion F comme la quantité de matière (particules) par 
seconde qui traverse l’unité d’aire d’une surface normale au mouvement de transfert 
étudié. F est aussi appelé la densité de courant de particules. En présence d’un 
gradient de concentration, on admet qu’il s’établit un flux de particules(un écoulement 
de particules) dans le sens descendant le gradient, et que ce flux est proportionnel au 
gradient correspondant: 

dC 

F = -D— = —DgradC, (1.2.14) 
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où D est appelé coefficient de diffusion ou diffusivité. 

Le signe négatif indique que le flux diffuse de la région ayant une forte concentration 
de particules à la moins forte. 

La seconde équation de Fick exprime qu’en tout point x la variation temporelle de 
la concentration C(x,t ) en fonction de sa variation spatiale au voisinage de ce point- 
par la relation: 


(1.2.15) 


de _ 

dt dx ’ 

cette équation décrit comment le changement dans la concentration dans un élément 
de volume est déterminé par le changement (variation) dans le flux entrant et le flux 
sortant du volume. 

En combinant les équations (1.2.14) et (1.2.15), on obtient: 


dC _ d D dC 
dt dx dx 


Si D est indépendant de la concentration, on obtient: 

dC _ d 2 C 
~dt ~ D ~d. 


Remarque 1.2.3 L’équation (1.2.14) est connue sous le nom de "Première Loi de 
Fick ". 


Remarque 1.2.4 La première équation de Fick (1.2.14) est utilisée en régime permanent, 
seulement, c ’est à dire lorsque le flux de particules ne dépend pas du temps. La deuxième 
loi de Fick exprime non plus un régime permanent de diffusion, mais un régime transitoire. 


1.2.3 Formes Quadratiques 

Soi la matrice réelle symétrique D= (aij)i<i,j< n ■ 


Définition 1.2.1 Le polynôme homogène du second degré relativement à n variables 

^1 ) ^ 2 1 i 'U*n 

n 

Y.aijUiUj. (1.2.16) 

i,j =1 

est appelé forme quadratique. 
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Désignons la matrice-colonne (u±,U 2 , .... u n ) par u et la forme quadratique par 


n 

A(u,u ) = aijUiUj , (1.2.17) 

i,j =i 


qu’on peut écrire 


A(u,u)=u t Du. (1.2.18) 

Théorème 1.2.2 Une forme quadratique est dite définie positive, i.e. 

A (u, u) > 0, u 0, 

si, et seulement si, tous les déterminants principaux successifs de sa m.atrice des coeffi¬ 
cients, sont positifs, i.e. 


det 1 > 0, det 2 > 0,..., det n > 0. (1.2.19) 

Preuve. Voir Gantmacher [13] page (308). ■ 


1.2.4 Quelques inégalités utiles 


Inégalité de Cauchy 

Soit a, b G M, e > 0, alors 


\ab\ < -a 2 + —b 2 
1 1 “ 2 2e 

Preuve. L’inégalité de Cauchy est un cas particulier de l’inégalité de Young suiv¬ 
ante: ■ 


Inégalité de Young 

Soit / une fonction continue et croissante sur [0, c], où c > 0. 
/ (0) = 0, a G [0, c] et b G [0, / (c)], alors 

pa pb 

ab< f(x)dx+ f~ l (x) dx 


où / 1 est la fonction inverse de /. 


( 1 . 2 . 20 ) 
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Preuve. On commence par l'expression 


g (a) — ab — / f(x)dx. (1.2.21) 

J o 

On prend b > 0 comme un paramètre. Puisque g' (a) = b — f (a) et la fonction / est- 
croissante, on a 

g' (a) > 0 pour 0 < a < / -1 (b) , 

g' (a) = 0 pour a = f~ l (b ), 

g' (a) < 0 pour a > / _1 (b). 

De cela, g (a) est une valeur maximale de la fonction g atteinte en a = f ( b ). 

Ainsi 

g (a) < ma xg (x) = g (/ _1 ( b )) . (1.2.22) 

Par intégration par parties on aura 


girHb)) = bf-'ib) 

rf-Hb) 




/ (x) dx 


(1.2.23) 


xf (x) dx. 

J o 

Si on prend y = f(x), l’égalité ci-dessus devient: 

9 (/ -1 ( 6 )) = / r 1 ( y)dy . 

J o 

En comparant (1.2.21), (1.2.23), et (1.2.22), on obtient (1.2.20) 

(Voir Mitrinovic, Pecaric et Fink [38]). ■ 

La fonction f(x) = x p ~ 1 avec p > 1 dans chaque intervalle [0, c] satisfait les condi¬ 
tions précédentes. On applique (1.2.20) utilisant —I— = 1 on obtient 

P Q 


a p b q 

Va, b G M + : ab <- 1 -. 

P 9 


(1.2.24) 

Si on remplace la fonction f(x) par ex p ~ l dans (1.2.20) où e > 0, alors on obtient 
l’inégalité de Young en e: 


Va, b G M + : ab < ea p + ^^b q . 

9 


(1.2.25) 


Ce qui donne 


Va, b G M + : ab < ea p + ( e)p ~ 1 b q . 
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Inégalité de Hôlder 

Soit p > 1 et q > 1 des nombres réels liés par la relation ^ ^ = 1 alors 

V (/,<?) e L»(n)xLi(n): 

1 1 

[ I f(x)g(x)\dx < ( [ \f(x)\ p dx] ( I \g (x)\ q dx) . (1.2.26) 

Jn \Jn J \Jn J 

1.2.5 Semi-Groupes et Générateurs Infinitésimaux 

Soit T(t ) une fonction définie sur [0, +oo[ à valeurs dans l’espace d’opérateurs linéaires 
bornés qui appliquent l’espace de Banach (E, ||.||) dans lui-même, T(t) : E —> E. 

Définition 1.2.2 On dit que T(t ) est un semi-groupe fortement continu ou de classe Cq 
si: 

i) T(0) = I. 

ii) T{t + s) —T(t) o T(s), pour tout t,s > 0, 

iii) T(t)u est continue comme fonction de t, pour chaque u fixée i.e. 

lim T(t)u = T(fo)u, pour tout to > 0 et tout u E E (fixée). 

t^to 

Théorème 1.2.3 Soit T(t ) un semi-groupe fortement, continu, alors il existe des con¬ 
stantes réelles u et C(uf) telles que 

||T(t)|| < C(uj) exput, pour tout t > 0 

Remarque 1.2.5 Dans le cas Cfu) = 1, lu = 0, le semi-groupe T(t ) est dit de contrac¬ 
tions. 

Définition 1.2.3 Considérons dans E le sous-espace vectoriel 

D(A ) = G E / fim — (T (A t) — I)u existe dans , 

et pour tout u G D(A), définissons l’opérateur A par 

Au = lim -r'- (T (A t) — /) u. 

At—>0 At 

Cet opérateur de domaine de définition D(A ) est appelé générateur infinit.ésim.al du semi- 
groupe T(t ). 
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Théorème 1.2.4 Si le semi-groupe T(t ) est de classe Cq, alors 

D(A ) est dense dans E. 

Preuve. Voir Pazy [41] page 5. ■ 

Proposition 1.2.1 Soit T(t ) un semi-groupe fortement continu sur E de générateur in¬ 
finitésimal A, alors on a 
i) Pour u G D(A); T(t)u G D(A) et on a 

————■u = AT(t)u = T(t)Au, pour tout t > 0. 

LLL 


U) Pour u E E 


t 

(T(s)uds G D(A) et A 
o 

Preuve. Voir Pazy [41], ■ 


T(s)uds 


= T{t)u — u. 


Définition 1.2.4 Soit A = [z G C : < argz < £ 2 , £1 < 0 < £ 2 ,} pour tout z G A 

soit T(z ) un opérateur linéaire dans un espace de Banach E. La famille T(z), z G A est 
un semi-groupe analytique (holom.orphe) sur A si 

i) z 1 —> T(z) est une application analytique sur A. 

ii) T( 0) = / et lim T(z)u = u pour tout u E E. 

z&K,z —>0 

ni) T(zi + z 2 ) = T(z 1 )T(z 2 ) pour zi, z 2 G A. 

Un semi-groupe est appelé analytique s’il est analytique sur un certain secteur A C C 
contenant le demi-axe réel positif 


Il est clair que la restriction d’un semi-groupe analytique à l’axe réel est un semi 
groupe continu. 

Considérons le problème d’évolution non homogène à valeur initiale 
f t u(t) = Au{t) + f(t) t t> 0 

■u(O) = Mo, 


(1.2.27) 


défini dans un espace de Banach E où / : [0, r) —» E et /l est le générateur infinitési 
mal d’un semi-groupe fortement continu T{t ) sur E. 
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Proposition 1.2.2 Si f G L x (0 ,t;E), alors pour tout uq G E, le problème (1.2.27) 
admet une seule solution donnée par 


u(t ) = T(t)u 0 + 



Preuve. Voir Pazy [41] page 106. 


— s)f(s)ds, 0 < t < t 
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Chapitre 2 


Introduction aux Systèmes de 
Réaction-Diffusion 

2.1 Introduction 

Les problèmes de Réaction-Diffusion s’écrivent sous une forme générale: 

du 

— - DAu = f(u), (2.1.1) 

où l’inconnue: u(x,t ) = (ui(x,t), ...,u m (x,t)) est un vecteur de variables, 
la réaction: f(x,t,u(x,t)) = t, u(x, t )),..., f m (x,t,u(x,t)) est généralement non 

linéaire, D est une matrice carrée m x m définie positive et diagonalisable appelée 
matrice de diffusion. Les termes de réaction sont le résultat de toute interaction entre 
les composantes de u, par exemple, en chimie u est un vecteur de concentrations chim¬ 
iques, et / représente l’effet des réactions chimiques sur ces concentrations. Le terme 
DAu représente les diffusions moléculaires à travers la frontière de la réaction. 

Les conditions aux bords seront choisies selon l’origine et la nature du problème 
étudié: s’il n’y a pas d’immigration des individus à travers la frontière du domaine 
D sur lequel le problème est posé, on choisit les conditions aux bords homogènes de 
Neumann. S’il n’y a pas d’individus sur la frontière, on prend les conditions aux bords 
homogènes de Dirichlet. L’inconnue (la solution qu’on cherche) est un vecteur dont 
les composantes sont généralement des fonctions positives: en chimie, par exemple, 
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c’est un vecteur de concentrations chimiques; en biochimie ou en métallurgie, c’est 
un vecteur de concentrations en nombres de molécules ou d’atomes respectivement; 
en dynamique des populations et en environnement, c’est un vecteur de densités de 
populations humaines, animales ou végétales... 

Les conditions initiales sont généralement positives; puisqu’il s’agit de concentra¬ 
tions, densités, charges électriques ... 


2.2 Modélisation 


Pour simplifier et exposer l’origine des équations de Réaction-Diffusion proposons nous 
à titre d’exemple l’étude d’une population avec densité u(x, t), vivant et déplaçant dans 
un conteneur. Pour décrire le mouvement, nous introduisons une autre quantité de 
charge; le flux de particules J(x, t ) G R n . A chaque emplacement x et à chaque instant 
t, le flux J(x,t ) est un vecteur qui pointe dans la direction générale de circulation à 
cet endroit. Sa grandeur |J(x, £)|, est proportionnelle à la quantité de particules qui 
circulent dans cette direction par unité de temps. Plus précisément, le flux J joue 
le même rôle que le flux de la chaleur, ou le flux de concentration pour un réacteur 
chimique. Considérons un volume d’essai Q de frontière T et nous équilibrons les flux 
entrants et sortants sur fl à travers T. Plus simplement: 

Changement de u dans Q = flux à travers r+ Changement relié à la naissance et la mort. 


Donc on obtient 

— f u(x,t)dV = — f J(x,t)dS+ f f(u(x,t))dV, 
dt Jçi J r Jn 

où f (u(x, t )) décrit le changement relié à la naissance et la mort, dV désigne l’intégration 
sur l’espace M n et dS représente l’intégration sur la surface M n_1 . 

Utilisant le théorème de divergence on aura 


d’où 


/ J(x,t)dS = / div J(x, t)dV, 

I r Jn 


—u — f (u) + div J ) dV = 0. 
>n \ dt 
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Cette dernière équation est satisfaite dans chaque volume de test Q. Ainsi (si la mesure 
dV n’est pas dégénéré) il en résulte que 

—u — f (u) + div J = 0. (2.2.2) 

CLL 

Ensuite, on a besoin d’une expression du flux en termes de distribution de population 
comme dans les réactions chimiques, ça nous conduit à utiliser la seconde loi de Fick 

J = —D'Vu. 

On suppose que le flux J est proportionnel au gradient négatif de la distribution de 
particules. Si nous combinons la loi d’équilibre (2.2.2) avec la seconde loi de Fick, nous 
obtenons l’équation de réaction-diffusion, 

—u = DAu + f (u). (2.2.3) 

dt 

Remarque 2.2.1 Dans le cas où f (u) — 0 , l’équation (2.2.3) devient l’équation de la 
chaleur. 


2.3 Exemples 

2.3.1 En bioécologie 

Le modèle décrivant la prédation des proies de densité spatiale u par des prédateurs 
de densité spatiale v est décrit par le système de réaction-diffusion 


du d 2 u du « 

a— — v — —h u(p — qv) = 0 


dt dx 2 


dx 


dv d 2 v du , . 

i dt ~ b w ~ “m + ” (r “ SM) = 0 


avec des conditions initiales vérifiant : 


dans (0, oo) x M, 


0 < 5 < Uq(x) < a 

0 < 7 < i; 0 (x) < f 


pour x G M, 
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où a, b > 0; p, r, q, s, u, /i, ô, a et [3 désignent des constantes positives et uo(x ) = 
■u(0,;r), Vo(x) = v(0,x ) sont des fonctions données. Les inconnues sont u = u(x,t ) et 
v = v(x, t). 

Physiquement, les constantes p et q sont les coefficients décrivant la croissance de 
la densité du prédateur, tandis que le coefficient .s décrit la décroissance de la densité 
de la proie due à la présence du prédateur. 

2.3.2 En électronique 

Prenons comme exemple la diffusion du phosphore dans le silicone. On dope une plaque 
de silicium par le phosphore (en la chauffant par exemple); le phénomène est décrit par 
les trois réactions: 


V + I 

< 0 > 

P + V 

E, 

P + I 

F. 


La première réaction décrit la neutralisation d’une lacune V (un lieu où se trouve 
un atome de silicium) par l’occupation de son lieu par un atome de Silicium interstitiel 
/ (c’est un atome de Silicium qui ne se trouve pas dans sa place habituelle). Dans 
la deuxième, un atome de Phosphore substitutionnel P prenant la place d’une lacune 
V, crée un complexe de type E, tandis que dans la troisième réaction, un atome de 
Phosphore substitutionnel P prenant la place d’un atome de Silicium interstitiel / crée 
un complexe de type F. Les trois réactions sons réversibles. 

Si on note par [V ], [E ], [/], [F] et [P] les concentrations en nombre de lacunes V, 
nombre de complexes du type E, nombre d’atomes de Silicium interstitiel /, nombre de 
complexes de type F et nombre d’atomes de Phosphore substitutionnel P , en appliquant 
la loi de conservation de la masse et en tenant compte de la diffusion de chaque type 
(par application de la loi de Fick J = —DVc, puis la loi de continuité VJ = — dc/dt ) 
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on aboutit au système de Réaction -Diffusion suivant: 

| [V] = d 1 A [V] - h [P] [V] + k 2 [E\ - k 0 ([R] [J] - [R e? ] [J^]) 

^ [P] = d x A [P] + k 2 [P] [R] - [E] 

[J] = d 3 A [J] - k 3 [P] [J] + h [P] - k 0 ([R] [J] - [V eq ] [J^]) , iGfi, t > 0 , 

| [P] = d 4 A [P] + /c 4 [P] [J] - k A [F] 

^ [P] = d 5 A [P] - h [P] [VI + k 2 [E] - k 3 [P] [J] + h [P] 

avec conditions aux bords bien choisies et conditions initiales positives. Voir Kouachi 

[ 26 ]. 

2.4 Existence globale des solutions 

Considérons comme exemple, le système de Réaction-Diffusion avec une matrice de 
diffusion triangulaire: 

u t = a Au — f (x, t, u, v), xÇ.Q,t> 0 , (2.4.4) 

v t = cAu + dAv — pf (x, t, u, v ), rrGf2,f>0 (2.4.5) 

supplémenté par les conditions aux limites et initiales suivantes: 
nu r)v 

^- = ^-=0, xedü ) t> 0 (2.4.6) 

au ou 

u (x, 0) = u 0 (x), v (. x , 0) = Vq {x) , iGÎÎ (2.4.7) 

où Q est un ouvert borné dans R N de frontière (El de classe C 1 . Supposons que: 

(Hl) Les constantes réelles a, d, c et p sont données par 

a > 0 , d > 0 , p > 0 et c G M. 

(H2) u 0 , v 0 G C (D) sont des fonctions non négatives. 

(H3) / : D x [ 0 , oo) xKxR —> M est une fonction continue mesurable localement 
lipshitzienne en u, v. 

(H4) / (. x , t, 0, r/) = 0 pour tout iGfi, t > 0, 7 / G M. 

(H5) / (. x , t , £, 77 ) > 0 pour tout x G O, t > 0, £ > 0, 77 G R. 
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(H6) il existe une fonction localement lipshitzienne ip dans R et deux constantes réelles 
k > 0 , a > 0 telles que 

0 < f (x, t, £, rj) < £;</?(£) e ari pour tout x G fi, t > 0, £ > 0, 77 > 0 

Le système (2.4.4)-(2.4.7) a été étudié extensivement dans les années récentes pour 
des cas variés, voir [2] [37] [17] [25] [23] [24] [3] [5]. 

Dans ce qui suit, on donne une réponse complète pour l’existence globale -en temps- 
des solutions de cette classe de systèmes et par une preuve unifiée, où nous utilisons la 
théorie des semi-groupes et la fonctionnelle de Lyapunov. 

Théorème 2.4.1 Si a ^ d, c > p (a — d) et v 0 — R 0; alors sous les hypothèses 

(H1)-(H6), il existe une solution locale unique nonnégative (u,v) pour le système (2-4-4)- 

(2.4. V- 


Preuve. Pour 1 < p < oc. les opérateurs A p et B p définis par 


D (A p ) = D(B p ) = lueW 2 ’ p (fl)- — \dfi =0 

A p u = a Au pour u G D (A p ), B p u = dAu pour u G D (B p ) 


génèrent deux semi-groupes analytiques sur l’espace L p (fi) notés S a ( t ) et S ( ] (t) re¬ 
spectivement. Par ailleurs, la restriction de — A p et — B p sur C (fl) est m-accretive. 
Soit E p l’opérateur défini par 

D(Ep) = \uGW 2 ’ p (fï)-,^\dfi =o|, 

E p u = cAu pour u G D (E p ), 


alors, l’opérateur matrice M p = 


A p 0 
E p Bp 


de domaine 


( _ 

Doc (Mp) = < u G W 2,p (fi) pour tout p > n, Au G C (fl) ; — \dfl = 0 

génère un semi-groupe analytique {S (£)} t>0 d’opérateurs linéaires bornés dans l’espace 
C (fl) x C (fl) (cf. [25] [41] et autres). 

Comme l’application, 


F (u (t ), v (t)) = (-/ (., t, u(t),v (t )), -pf (.,t,u (t ), v (t))) 
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est continue mesurable localement Lipschitzienne en u. v dans l’espace C (fl) x C (fl), 
alors le problème (2.4.4)-(2.4.7) admet une solution classique locale unique dans l'intervalle 
(0, T) . Par ailleurs, il est bien connu que l’implication suivante est vraie 

(Vf G ( 0 , T max ), ||u (t)^ + ||u (f)^ < Const .) => T max = +oo, 

où T max est le temps maximal d’existence 

Par le principe de maximum (cf. [44]), il est clair que d’après la condition (H3) on a 
0 < u < KL pour tout x G fl, t > 0 (2.4.8) 


• Si c = 0, et comme u > 0, alors / (x,t,u,v) > 0 d’après (H5) et par (H6) et le 
théorème de comparaison (cf. [44]), on obtient v > 0. 

• Si c > p (a — d) et vo — > 0, en posant w = v — , alors d’après (2.4.4)- 

(2.4.7), on aura 


w t 

dw 

du 

w 0 


dAiv + 



/ (x, t, u, v) sur fl x ( 0 , oo), 


0 sur <9fl x ( 0 , oo ), 
c 

V 0 - zUo- 

a — d 


(2.4.9) 


Par le théorème de comparaison (cf. [44]), on obtient que ru > 0, et par conséquent 


v > ——-u > 0 
a — d 


Théorème 2.4.2 Sous les hypothèses du théorème 2-4-1, la solution du système (2-4-4)- 
(2-4-7) existe globalement. 

Pour démontrer ce théorème, commençons par le lemme suivant, 
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Lemme 2.4.1 Supposons a d . Soit (u,v) la solution du système (2.4-4)-(2-4-V donnée 
par le Théorème 2-4-1; alors, pour toute constante 7 > 0, il existe une constante (3 = 
f3 ( 7 ) > 0 telle que la fonctionnelle 


1 1 —> L (u (t), v (t)) = / e (u+ ÏÏ 2+ ^’dx 


est noncroissante sur (0, T max ). 

Noter que L est une fonctionnelle de Lyapuvov. 

Preuve. Notons w = (u + /3) 2 + 'yv. Les équations (2.4.4)-(2.4.5) donnent pour 
tout t G (0, T max ) 


m = 


2 (jï) (u +/3)+ 7 (§?) ) eiu+p)2+lvdx 

(2 (u + (3) (aAu — f) + 7 (cAu + dAu — pf)) +TU do; 


donc 


= 2 a (u + (3) A ue w dx — 2 / (u + /3) fe w dx + 7 c / A ue w dx 

Jfi Jn Jn 

+7d f A ve w dx — 'ypf fe w dx 

J Çl J ÇL 


L\t) = 2a (u + /3) Aue w dx + 7 c / A ue w dx 


+7 d / A ve w dx + / [— 7 /? — 2 (w + /?)] fe w dx 

Jn J fl 

= U + ^2 + ^3 + I 4 . 


(2.4.10) 


Utilisant la formule de Green, on aura 


U = 2 a (u + f3) A ue w dx = —2a / VuV [(u + /5) e™] dx 


J fi 

= —2a / Vu 

J fi 


(Vu) e™ + (Ve (u+/3)2+7,;N ) (u + /3) 


da; 
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= -2a Vu, {(Vu) + (u + /3) [2 (u + (3) (Vu) + 7 Vu]} e w dx 

Jn 

n 


- 2 »e / 

7=1 


du f du 


— Jndxj [dxj 


+ (u + /3) 


^ du dv 

2 {u+f 3 ) dï j + ' ï & Cj 


e w dx 


= -2a [ {| Vu | 2 + {u + (3) [2 (u + (3) | Vu | 2 + 7 VuVu] } e w dx 
Jn 

= —2a f | Vu | 2 e™<7r — 4a f (u + /3 ) 2 | Vu | 2 e w dx 


-2 a 7 / (u + /3) ('Vu) ('Vu) e™c7r 

Jn 


(2.4.11) 


et 


/ 2 = 7 c / Aue w dx = - 7 c / (Vu) Ve™c7r = - 7 c / (Vu) Ve ( “ +/ 3 ) 2 + 7 î; cfa; 


— 7 c / (Vu) [2 (u + /3) (Vu) + 7 Vu] e™<ir 


v-^ /‘ <9u 

' 7C ^- / , 

J=1 ax j L 


_ <9u dv 

2 ( “ + 0) + 7 a7 


e w dx 


—2 7 c / (u + /3) |Vuj" e w dx — 7 2 c / (Vu) (Vu) e"’<7r 

Jn Jn 


(2.4.12) 


et 



7.3 = 7 d / Auewdx = - 7 c/ / (Vu) Ve™du = - 7 d / (Vu) Ve (u+/3) + 7 ^a; 

J r2 J fi J 

u + /3 ) 2 + 7 u) +ry c7c 

e w dx 



_ s <9u dv 

2{u + p) W j +1 â^ u 


= — 270 ? / (u + /3) (Vu) (Vu) e™c7r - 7 2 d / |Vu| e™<7r. 
7 Q ./n 


(2.4.13) 


Comme f3 < u+(3 < ||uo||+/3, alors par l’inégalité de Young-e (|aé| < ^a 2 + e& 2 ) , on 
obtient 


Ii =—2a / |Vu| e™<ir — 4a / (u +/3) |Vu| e w dx — 2ar/ / (u + /5) (Vu) (Vu) e™<7r 

J Çl J fl J f2 
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< —2a / | Vu | 2 e w dx — 4a / (3 2 | Vu,\ 2 e w dx + 2a 7 / (u + /?) |(Vu) (Vv)| e w dx 

J fi J fi J fi 

< -2a [ | V| 2 e w dx — 4a [ (3 2 | Vu\ 2 e w dx + 2 a 7 (||u 0 || + P) I |(Vu) (Vu)| e w dx 

J Çl J Q J Çl 

< —2a f | V «| 2 e w dx — 4a f (3 2 | Vu) 2 e w dx + 2a 7 (||« 0 || + /3) f( (Vw )| 2 + e |(V'o )| 2 ) e 


h < 


Ainsi 

—2a — 4 a/3 2 
+2a7(||«o|| +P) j~ e 


4e 


\Vu\ 2 e w dx + 2a^(\\u 0 \\ + (3) e / \(Vv)\ 2 e w dx (2.4.14) 


et 


/ 2 = — 27 c / (u + (3) |Vu| ’ e w dx — 7 2 c / (Vu) (Vu) e w dx 

J fi J fi 

< 2-y |c| f (u + f3) \V u\ 2 e w dx + j 2 \c\ f \(Vu) (Vv)\e w dx 

J fi J fi 

< 2 7 |c| [ (||«o|| + /?) |Va| 2 + 7 2 |c| [ (\- |(V«)| 2 + e |(Vv)| 2 ) e w dx. 


4e 


Ainsi 


V 


h < ( 2 7 |c| (||tt 0 || + /3) + 7 2 \c\ — ) / |Vu|“ e w dx + 7 2 |c| e / |(Vt>)| e w dx (2.4.15) 


et 


/ 3 = -2 7 d / (a + /3) (Vu) (Vu) e w dx - 7 2 d / | Vv| 2 

J fi J fi 

< 2 7 d [ (u + (3) \(Vu)(Vv)\e w dx--/ 2 d [ \Vv\ 2 e w dx 

J fi J fi 

12 


< 2: 7 d / (K|| + /?) ( ^|(Vîi)r + e|(ViOn e w dx- 1 2 d / |V< e w da; 


alors 


J 3 < 2 7 <i (||m 0 || +/3) — / | Vu | 2 + [2 7 d (||« 0 || + /?) e — 7 2 rf] /"|Vu| 2 e w rfa:. (2.4.16) 

Il est clair que /4 = f Q [— 7 p — 2 (-u + /?)] fe w dx est négative. 

Donc de (2.4.14)-(2.4.16) et (2.4.10) on trouve 


m < 


—2a — 4a/3 2 + 2a 7 (||'Uo|| + /3) — 


|Va| 2 é w dx + 2a 7 (||'U 0 || + /?) e / |(Vv)| 2 e w dx 


2 7 |c| (||m 0 || +/3) + 7 2 |c| — / |Vu,| 2 + 7 2 |c| e / |(Vv)| 2 e w dx 

4e J Jfi Jfi 

+2 7 d (||m 0 || + /?) t 1 [ \Vu\ 2 e w dx+ [2 7 d(||a 0 || +/3)e- 7 2 d] f | Vv| 2 
4e Jn Jo 


dx. 
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ou encore 


L'(t) < A / |Vîx| : e w dx + B / \{Vv)\ 2 e w dx 

J Q. Jn 


(2.4.17) 


ou 


A = 


-2a — 4a/3 2 + 2a7 (H^oll + /3) —— 


4e 


+ 


2 7|c| (Kll + P) +7 2 |c| ^ 


+ 2 7 C KKII + k) ^ 


(2.4.18) 

B = 2 a 7 (Kll + /5) e + 7 2 kl e + [ 27 c? ( 11 xto 11 + P) c — 7 2( ^] • (2.4.19) 

Estimons e>0, 3 > 0 et 7 > 0 pour avoir /I < 0 et B < 0. Pour cet object, /I < 0 


si 


e > 1 2o 7 (Hti.ll + /5) + 7 2 kl + 27 c? (Kll + /3) 


8 a + 2a/7 -7 kl (Kll + Æ) 


Aussi, 7? < 0, si 


e < 


7<7 


2 (a + rf) (Kll + P) + 7 l c l 

Choisissons P pour obtenir 

1 207 (Kll + /5) + 7 2 kl + K (INI + Æ) 


8 a + 2a/3 2 — 7 |c| (|K|| + /?) 


< 


7 <7 


2 (a + h) (11i/o11 + P) + 7 l c l 


ou 


4 (4a<7 — (a + h) 2 ) /3 2 — 4 (2 (a + d) ||u 0 || + |c| ( a + 3rf) 7 ) P 


— |c | 2 7 2 — 4 |c| (a + 3rf) Kll 7 — 4 (a + <7)"“ 


A 8a<7 


> 0. 


Cette inégalité est réalisée si on choisit P de manière que 


( 2 (a + h ) 2 \\u 0 \\ 
y + |c| (a + 3 d) 7 


/?< 


/ 


16 


-16 (4a<7 — (a + d) 2 ) 


2 (a + d ) 2 KH 
+ |c| (« + 3<7) 7 

c| 2 7 2 + 4 |c| ( a + 3d) Kll 7 

4-4 (a 4- dy KH "b 


8 (4a<7 — (a 4 - (7)") 


\ 


1/2 


/ 
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Donc, si nous prenons 


1 207 (||m 0 || + fd) + 7 2 H + 27 d (K|| + fi ) < < _ jd _ 

8 a + 2a/3 2 — 7 |c| (||«o|| + f3) 2 (a + d) (||rco|| + /?) + 7 |c 


/ 2 (a + d ) 2 KH 
y + |c| (a + 3c?) 7 


P< 


( 


2 (a + d) 2 ||« 0 || 

+ |c| (a + 3 d) 7 

(iad-(a + df)\ HV + 4|c|(« + 8d)|N| 7 

+4 (a + d ) 2 ||m 0 || 2 + 8 ad 


\ 


1/2 


2 (4acl — (a + cl) ) 


on obtient que L'(t) < 0, pour tout t E (0, T max ) ■ 

Preuve du Théorème 2.4.2. 

Comme est une fonction continue et comme u est bornée; alors, à partir de la 
condition (H 6 ), il existe une constante positive M > 0 telle que 

0 <f(x, t, u , v) < Me av . (2.4.20) 


Choisissons 7 = na dans le Lemme 2.4.1; alors il existe une constante [3 > 0 telle que 
pour la fonctionnelle 

L 0 (it ) = [ e {u+ P )2+nav dx (2.4.21) 

Jn 

on aura 

I/o (t) < 0, pour tout t E (0, T max ) . (2.4.22) 

Comme u > 0, alors de (2.4.21)-(2.4.22) on obtient 

[ e nav ( t '>d x < meas (O) e (^°^ +/3 ) +n °'lbolloo pour tout t E (0, T max ). (2.4.23) 

Jn 

Il en résulte que e av E L°° ([0, T max ) ; L n (D)) et 


sup ||e CTî; ^|| < [meas (Q)] 1//n e n ^ U0 ^° o+l3 ') +a ^ vo ^°° 

0 <C<T max 


(2.4.24) 


et par suite, de (2.4.20) et (2.4.24), on aura /(.,., u, v) E L°° ([0, T max ) ; L n (D)) et 


sup ||/ u, u)|| n < M [meas (0)] 1/n e«0' Uo ^°° +/3 ) +<T ^ 0 ^°° (2.4.25) 

0 <C<T max 
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où ||.|| n est la norme de l’espace de Lebesgue L n (fl). 

De (2.4.9), on obtient 

v ( t) = S d ( t ) ( v 0 - C —J u 0 ) H- —jU (t) 

\ a — a J a — d 

+ (^~P+J o S d (t - s) f (.,s,u(s) ,v(s))ds. (2.4.26) 

Finalement, pour a>0et0<T — a < T max , nous avons à partir de (2.4.9) que pour 
tout t E (0, T) 

v(t + a ) = S d (a) v ( t ) H- —u (t + a) - —.Sd (a) u (t) 

a — a a — a 

+ p + J S d (t — s) f (^.t + a — s,u(t + a — s) ,v (t + a — s)) ds. 

(2.4.27) 


Utilisant l’estimation (1.1) de [16] avec s = |, il en résulte qu’il existe une constante 
c (n) > 0 telle que si 0 < a < 1 , alors 

l|s J («)W0IL<^5ll»’WL, 

Il S d (t - s) f (,,t + a - s,u(t + a - s) ,v (t + a - s))^ 

C ( 77 ') 

< -jzÿk 11 / (-,t + a - s, u (t + a - s) , v (t + a - s))|| n . 

Comme v n < e nav , nous obtenons par (2.4.24) et (2.4.25) que 

\\Sd (et) v (i)|| < [meas( 0 )] 1 /n eHlNII-+^) 2 +^olU ) (2.4.28) 

a z/ d na 

\\S d (t - s) f (.,tS- a - s,u(t + a - s) ,v (t + a - s))^ 

< ^M[meas(0)] li/n e™( ||uo|l - +/3 ) 2+CT||î ’ 011 -. (2.4.29) 

s l/ A 

En outre 

(2.4.30) 


où M d et c j sont telles que \\S d (t)|| < M d e üja , pour tout t > 0. 
Via (2.4.28), (2.4.29) et (2.4.30) et par (2.4.27), on obtient 


\v (t + a) 


< 


a — d 

+c (■ n ) 


IKIL (1 + M d e“ 


n\ 


a 


2/3 


TUT 


+ 3 M\fâ f-p + 


V 


a — d 


[meas (f ï)] l ^ n e n 0l'“°ll oo +/3) + cr ll' L ’oll c 


(2.4.31) 



d’après (2.4.31) on aura 


P Wlloo < c ( a > IMLi'IMU pour tout t G [0,T max ). (2.4.32) 

Par conséquent, à partir de (2.4.8) et (2.4.32), on établit l’existence globale de la 
solution (u,v); c’est-à-dire T max = +oo. ■ 
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Chapitre 3 


Systèmes de Réaction-Diffusion de 
type Gierer-Meinhardt 

3.1 Introduction 

La formation des modèles dans le développement d’une organisation plus élevée hors 
d’un œuf simple fertilisé ou d’un tissu végétal est l’aspect le plus fascinant de la biologie. 
Une question centrale est pourquoi les palmes se disposent d’une manière organisée sur 
le stipe d’une plante? Cependant, un regard sur la nature inorganique indique que la 
formation de modèles (En anglais: Pattern formation) n’est pas particulière aux objets 
vivants. La formation de modèles est également la règle dans le monde non vivant; 
comme la formation de galaxies, nuages, chute de pluie, foudres, cristaux, érosion, 
toutes témoignent de la génération de structures commandées. 

Il est instructif pour rechercher des principes communs dans la génération de ces 
structures; si une petite déviation d’une distribution homogène a une rétroaction posi¬ 
tive forte sur elle-même, la déviation augmentera. Par exemple, l’érosion procède plus 
rapidement à l’endroit ayant déjà subi des dommages initiaux aléatoires, et les fleuves 
contournés sont formés en dépit de la réparation homogène de la pluie presque tout le 
pays. Une grande dune de sable peut résulter d’une pierre dans un désert qui produit 
un abri de vent et peut accélérer ainsi localement le dépôt du sable; ce dépôt augmente 
l’avantage de la probabilité de dépôt de sable, et ainsi de suite. 
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En Biologie; Alfred Gierer et Hans Meinhardt ont formalisé cette observation et ont 
proposé un modèle moléculaire plausible pour modéliser la formation, composé de deux 
équations aux dérivées partielles formant un Système de Réaction-Diffusion à réaction 
fractionnaire de type activateur-inhibiteur: 


d a a 2 d 2 a 

Wt = P J - + a + Pa ^ ^ ^ £ ^ ( > „ 

- =p a*-, h h + D h — 2 


(3.1.1) 


Dans ce système, a est une substance autocatalytique à courte portée, (i.e.) activateur 

Q a 

et h est son antagoniste à longue portée, (i.e.) inhibiteur, — décrit le changement 

de la concentration de l’activateur a par rapport au temps. Le premier terme du 
a 2 

membre de droite (p—) décrit le taux de production qui dépend non linéairement de 

la concentration de l’activateur et de l’inhibiteur. Le nombre de molécules délabrées 

par rapport au temps est proportionnel au taux de décomposition (p a ) et au nombre 

de molécules actuelles (a) (comme dans une ville le nombre de personnes décédées 

dépend du nombre d’habitants). L’échange des molécules est modélisé par la diffusion 
d 2 a 

(D a -—), (il y a d’autres modes possibles de propagation). La deuxième équation décrit 

(J ■' 

en termes analogues le changement de la concentration de l’inhibiteur. p a est un petit 
taux de production d’activateur qui est exigé pour lancer l’autocatalyse d’activateur à 
une concentration très basse (par exemple, dans le cas de la régénération). 

Le modèle décrit la concentration d’une substance autocatalytique à courte portée; 
l’activateur, qui règle la production de son antagoniste à longue portée, l’inhibiteur 
(Gierer et Meinhardt, 1972, Gierer, 1981, Meinhardt, 1982). C’est certainement un 
modèle minimal, mais il fournit un pont théorique entre les observations d’une part et 
la déduction des mécanismes moléculaire-génétiques fondamentaux d’autre part. 

La possibilité de construire de modèles par la réaction de deux substances qui dif¬ 
fusent avec différents taux a été découverte par Turing (1952) [48]. Gierer et Meinhardt 
[15] ont prouvé que même si cette condition est satisfaite, seulement une classe très 
spéciale des réactions peut former des modèles dans le cas de l’autocatalyse locale 
et l’inhibition de long-portée. Des différentes réalisations sont possibles : L’inhibition 
peut résulter d’un épuisement où l’autocatalyse peut être basée sur une inhibition d’une 
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inhibition. L’équation de Turing satisfait cette condition. Au terme de mécanisme de 
Gierer-Meinhardt, l’inhibition résulte d’un déplacement d’activateur par la diffusion 
rapide de la substance. 


3.2 Exemple de modélisation 


Comment la coquille de mer s’embellit avec de telles couleurs? 

Une cellule productrice-colorante infecte une voisine et après une partie de retard 
cette cellule également devient coloré et déclenche alternativement la prochaine voisine, 
et ainsi de suite jusqu’à la génération des lignes obliques. 

Les modèles de la production de colorants sont plus communs avec les modèles 
des formations. La situation générale presque dans tous les modèles qui forment des 
systèmes est que les groupes des cellules se développent différemment de leurs voisines. 

Pour la modélisation, supposons que la coloration dans la coquille est sous la 
commande d’une molécule; l’activateur a. La réaction inhibitrice peut résulter de 
l’épuisement d’un substrat, s. 

La formulation mathématique de ces interactions (activateur-substrat) est donnée par 
le système de réaction-diffusion de type Gierer-Meinhardt (compliqué) suivant (voir 
Meinhardt [31]): 


d a 9 * d 2 a 

— = psa z — pa + D a —— 


d s 
dt 


d t 


= a — psa 


dx 2 l 
d 2 s 

- z/s + D s 

ox z 


où 


1 + ka 2 


+ Po- 


L’activateur, a et sa molécule précurseuse, le substrat s, diffusent d’une cellule à une 
autre avec les taux D a et D s respectivement et se délabrent avec les taux p et v. 
Le substrat s produit avec le taux constant a. La constante k est un paramètre de 
contrôle, p représente la densité de source et p 0 est une valeur initiale de production. 
Ces équations nous permettent de calculer le changement de concentration de a et de 
s pendant un intervalle de temps. 
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3.3 Etude de l’existence globale des solutions de 
système de Gierer-Meinhardt d’ordre deux 


Le système de Gierer-Meinhardt (3.1.1) d’ordre deux est l’origine des systèmes formés 
d’un activateur (a) et inhibiteur (h) , l’existence globale des solutions avec conditions 
aux bord de Neumann et conditions initiales a ( x , 0) > 0 et h (x, 0) > 0 a été prouvée 
par Roth [44] en 1984 pour N < 3 , le même résultat est obtenu par Wu et Li [49] en 
1990, à condition que a, h x et p a soient convenablement petits. 

D’une manière générale considérons le système de Réaction-Diffusion qui généralise 
le système de Gierer-Meinhardt (3.1.1) (voir [15]) 


dui . , 

— = OiAui - fiu\ + b 1 -^ +a, 
(JL 

du 2 


= a 2 Au 2 - vu 2 + b 2 ^~, 
dt u s 2 


iGfi,t> 0, 
x G O, t > 0, 


(3.3.2) 


avec conditions aux bords de Neumann 


du i du 9 

—- = 0 , —- 
d r] ’ d r] 


0, 


et conditions initiales 


x G T, t > 0 , 


ui(a:,0) = ^i{x) > 0, iefl, 

«2(^,0) = i/j 2 (x) > 0 , x G D, 


(3.3.3) 


(3.3.4) 

(3.3.5) 


où: 

C est un domaine borné de frontière T régulière. 

Les constantes de diffusions: a± et a 2 sont supposées positives. 

/i, u, a, bi et b- 2 sont des constantes positives, 
p, g, r et s: des indices non négatifs avec p > 1. 

0 < 61 < b- 2 . 

En 1987, K. Masuda et K. Takahashi [36] ont prouvé l’existence globale des solutions 
du système (3.3.2) sous la condition . Après quelques années, L. Mingde, C. 

Shaohua et Q. Yuchun [35] ont trouvé en 1995 de meilleurs résultats, ils ont démontré 
l’existence globale sous la condition r > p — 1 et rq > (p — 1) (s +1) pour des conditions 
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initiales arbitraires, ils ont prouvé aussi l’explosion en temps fini des solutions dans 
le cas où r > p — 1 avec rq < (p — l)(s + 1 ) ou r > p — 1 pour certaines valeurs 
initiales, presque les mêmes résultats ont été démontré par Ph. Souplet et P. Quitner 
[43], Ils ont établi l’existence globale des solutions sous la condition V < min (iTÏ,l) 
et l’explosion sous > min l) , fi 1 . 

Remarque 3.3.1 1) Les auteurs de [35] et [43] ont prouvé l’explosion en temps fini pour 
les systèmes EDO associés. 

2) Dans les deux derniers travaux, les auteurs ont supposé que bi — b 2 — 1, on va 
analyser ces méthodes mais avec une hypothèse plus générale; 0 < b\ < b 2 . 

3) Pour plus de détail sur la discussion d’existence locale des solutions, voir [18], [45], 
[43] et autres. 

3.3.1 Principe de Maximum 

Les résultats essentiels sont les suivants 

I) Soit («!(£,.), u 2 (t,.)) la solution de (3.3.2). Alors pour tout (t,x) dans (0, T max ) x 
fl, on a 

Ui(t,x) > e~ blt min (fifix)) > 0, (3.3.6) 

o 

u 2 (t,x) > e~ b2t min > 0. (3.3.7) 

o 

Ces bornes nous donne que si T max < 00 , alors fini sup ||rq = 00 

t ^max 

II) Il existe une constante c > 0 telle que 

Ml, «2 > c, x G flj 0 < t < T max . (3.3.8) 

Comme (ui) t — a^Aui > 0 dans {«i < o/fi}, on a U\ > := min {(J/p, min^ fifi) 

dans fl x [0, oo[. Alors u 2 satisfait (u 2 ) t — a 2 Au 2 > 0 dans |«2 < (S[/n) 1/( ~ A+1 ^ | et 
u 2 > S 2 := min ( \8\/v) l/ <s+1 ^ ,mhijy^ 2 j dans fl x [0,oo[. Prenons c = max(^i,5 2 ) • 
voir [43]. 
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Lemme 3.3.1 Soit p, q,r et s des constantes telles que, r > p — 1 et rq > (p — l)(s + 1). 
Pour tous h, a, f3 > 0, il existe c = c(h, a, (3) > 0 et 9 = 9 (a) G (0,1), telles que, 


p- 1+0 V r+a /v a \ u 


Preuve, (voir [43] page 298) 


u > 0 , v > h. 


Lemme 3.3.2 Supposons que u(x ) > 0 et v(x) > 0. Alors pour tout ensemble d’indices 
p 0 , q 0 , a, (3 et A satisfaisant À < p 0 < a (non nécessairement positifs) et pour toute 
constante c > 0 , on a 


— dn <C —jltt + c -(P 0 -A)/(a- P0 ) / “ 


6 = [q 0 (a - À) - /3(p 0 - A)] (a - p 0 ) 1 . 

Preuve. Par une application de l’inégalité de Young on achève la preuve.(voir 
[35]). ■ 

Lemme 3.3.3 Soient u±, « 2 , solutions du problème (3.3.2). Alors on a: 


— < max (^ Bi. 


kl - wïïmJo» 

il - 


Preuve. 1) On a 


— [ —rdQ = —k [ —r—r (a\Aui — uu\ + 6 iK + o j d£l, où k est un entier positif, 
dtj uf J ul +l V U 2 J 

n a 

/ I Ç 1 Ç 11 F —^^ Ç 1 

... AuidQ + kn / —rdPl — kbi / -A—- — dit — ka / ... dft. 

«i J uf i m] J u f +1 
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Appliquons la formule de Green pour le terme 




Auidfl, (en tenant compte des 


conditions aux bords de Neumann), et puis utilisons l’inégalité de Young pour le terme 


kfi J —dn pour obtenir 
o 1 


4- [< ka~ k n k+1 101 . 
dtj u\ ~ 11 

fi 

Intégrant les deux membres de cette inégalité de 0 à t, on aura 


4 dO< [~4dn + ka~ k / +1 \n\t. 


Uî 


r 


En élevant les deux membres à la puissance \ et en faisant tendre k vers +oc, on trouve 


1 

Ul 


^ (k 

< max( — 
o 




= B X . 


2) D’une façon similaire on a 


d r î 

j t h/n=- k ( k +i) a 2J 


1 |V« 2 | 2 <iSÎ + kv [igdU - kb 2 [ -+i dÇl ' 

u 2 


Un 


En utilisant le Lemme 3.3.2 pour / —rcZO, on obtient 

J «2 

n 




u 


V J Un 


V 


— k/s+1 


kr/s+1 

«1 


dn, 


ce qui donne 


1 k ( 1 \ «+1 1 1 k / 1 \ s +1 1 

-dn<k(b 2 )^\^-j j —— dü<k{b 2 ) N 

fi 1 


i 

Ul 


kr/s+1 


Intégrons de 0 à t, élevons à la puissance \ et faisons tendre k vers +oo, pour avoir 

1 


1 

u 2 


< max | | — 
v 


bn \ S+1 — 

2 \ DS + 1 

n l i 


lp 2 


= Bo. 
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Théorème 3.3.1 Supposons que r > p — 1 et rq > (s + 1) (p — 1) , alors, pour toute 
V’i, ^ 2 » V’ï 1 »^ 1 e T°°(0)j le problème (3.3.2) admet une solution globale et bornée dans 
l’intervalle de temps (0, oo). 


w\t) 


Preuve. Pour tout n > 2, soit rn > 0 suffisamment petit, tel que 
mn 2 (ai + a 2 ) 2 


4 (m + l)a 2 


< apn(n — 1), et n/i — 3 mu > 0. 


Soit «,(*) =ju > ~ 2 -dn. Démontrons d’abords que «**) est bornée pour tout t > 0. 
n 

Dérivant w(t), on obtient 

f U l~ l ( A U 1 \ JO 

i / — 1 a iA«i ~ A^i + Oi-g + cr 1 dD - 

J u 2 V U 2 J 


= n 


U ( . , M \ 

m 1 Dïï+ï t « 2 Au 2 - ^«2 + » 2 — ) aD 


fu n ~ 2 f u n ~ l 

n(n — l)oi / — 11 — |Vwi| dû + mn{a\ + a 2 ) / ,^ +1 VuiVu 2 dO — fimw{t ) 


ut 


u, 


+ra&i / ?il m+q + ncr / dO - m(m + l)a 2 / ^+ 2 | Vu 2 | 2 dD + umw(t) 


n 

n—l+p 


U 


,n—l 


U V 


U 


m-\-q 


Un 


U 


-rribÿ 


u 


n+r 


m+s+1 

a 2 


dn. 


Utilisant l’inégalité de Cauchy on aura 


mn(ai + a 2 ) / . VuiVu 2 dO 


n—1 




m +1 


< m(m + l)a 2 / | Vu 2 | 2 d(2 


u 


?n+2 


mn 2 (ai + a 2 ) 2 f u\ 


n —2 


4(m + l)a 2 


ut 


| Vui| 2 dfh 


Donc on a 


n—1+p 


n+r 


n—1 


w'(t) < nbi / —dD - mb 2 / - (/m - vm)w(t ) + rwx / dfi. 


u 


u 


ut 


n n 

Après, choisissons e G (0, n) tel que 


t= nh-(» + l)(p-l)] +ii ,-,-l +m a 


e(r + 1 — p) 


r + 1 — p 
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appliquant le Lemme 3.3.2, en prenant a — r + n, /3 — s + l + m et \ — n — e, on trouve 

, f u r 1+ï \ 0 ^ , f < +r 


0 = [(m + g) (r + e) — (s + 1 + m) {p — 1 + e)] (r + 1 — p) 1 
= [rg — (s + 1) (p — 1) + e (g — s — 1)] (r + 1 — p ) -1 + m, 

Çl = n6i(m6 2 /n6 1 )- (p - 1+e)/(r+1 - ï,) . 

Ensuite, appliquant le Lemme 3.3.2 en prenant: a — n, /3 — m e t A = 0, on aura 

fn n ~ e C v n /‘I 


t / e)/e 

ç 2 = Çi {mu/ çi) . 


D’une manière analogue, on obtient 


Enfin, on a 


fu n ~ f u n f 1 

na / ——dil < mu / — dU + ç 3 / —c/O, 

I „ .777. - I „ . 777, ° / „ . 777, ' 


ç 3 = n [crmz// (nu)] (r! 


«/(£) < Ç 2 [~üdü + ç 3 [-\rdQ - (■np - 3 mu) w(t), 


multipliant par exp [(np — 3 mu) t] et appliquant le Lemme 3.3.3, on obtient 
— (w(t) exp [(np — 3 mu) t]) < (ç 2 L ?2 + Ç 3 5™) |0| exp [(np — 3mn) t]. 

LLL 

Intégrant de 0 à t, on aura 


w{t) = / —^dQ < (ç 2 i ?2 + Ç 3 B™) |0| exp [(np - 3mn) *] + w(0) := B 3 . (3.3.9) 

./ Ur. 


38 



Maintenant, démontrons que U\ et 112 sont bornées en utilisant l’estimation (3.3.9). 
Il est bien connu que le problème (3.3.2) a une paire de solutions locales uniques 
Ui, u 2 G C 2,1 (II x (0,2t)) , pour certain r > 0. Supposons Au\ = dAu x — /mi, alors 
A est le générateur infinitésimal d’un semi-groupe holomorphe { e tA |0 < t < 00 } dans 
L a (H), de domaine 

D(A) = ew 2 > a = o\, 

où a = IV (N + 1) /2 et N est la dimension de M. N . Par ailleurs il existe des constantes 
ç > 1 et 5 > 0, telles que 


A 2 e 


t.A 


< Ç ( 1 + t 


-St 


t > 0 , 


(3.3.10) 


où 
du 1 
dr] 


dénote la norme de L a (H). Ainsi la première équation du système (3.3.2) avec 
0 et Ui (x, 0 ) = ^ 1(^)5 est équivalente à l’équation intégrale 


t 

ui ( t ) = ê A t\) 1 + Je^~ s ^ A [b\Ui (s) u^ q (s) + a ] ds. 
0 

Utilisant le théorème 1.6.1 cité au [18] et (3.3.10), on obtient 


nraxM 1 (.,t) < A^u\ (t) 


< ç (1 + H) e~ 5t ll^ll + ç f e-W-^ (t - s) 


, u i 

kl ~ g S" CT 

u\ 


(3.3.11) 

ds. 


Supposons n > ap et nq > pm. Appliquant le Lemme 3.3.2, on trouve 


U\ 


Uo 


= f%dTl< [A d{l + IA da , 

J u\ J U™ J u\ 


(3.3.12) 


où 7 = (aqn — apm ) (n — ap) 1 > 0. Alors les estimations (3.3.9), le Lemme 3.3.3, et 
les équations (3.3.11), (3.3.12) impliquent 


rriaxîi! (., t) < M 


1 + t 


t > T. 


(3.3.13) 


Il résulte de (3.3.13) et du principe du maximum que la solution est bornée. ■ 
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Chapitre 4 


Existence globale de solutions pour 
certains systèmes de R-D de type 
Gierer-Meinhardt 


4.1 Existence globale et Explosion en temps fini de 
solutions à trois composantes de systèmes de 
R-D de type Gierer-Meinhardt 


Le développement d’une plante, en biologie végétale, passe par une chaine d’étapes de 
formation des modèles, et parmi ces étapes, le phénomène de disposition des palmes 
sur le stipe (la tige) de la plante qu’on appelle en terme biologique: La Phyllotaxie. 
Meinhardt, Koch et Bernasconi [33] ont fait des études de modélisation dans ce do¬ 
maine et ont proposé un modèle mathématique qui décrit un mécanisme de réaction 
fractionnaire entre un Activateur a et deux inhibiteurs h et s: 


d a d 2 a 

dt = Da d.X* - /ia ° + h(s + k a ) 
dh n Ô 2 h 2 

dt = Dh w~ ,ihh+a * 

ds ers 

dt = Ds d^ ~ ^ + 


+ (?a, 


(4.1.1) 
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Plus généralement, considérons le système de réaction-diffusion à trois équations 

(du u pi 

— = ai A u - p 3 u + + a, x G il, t > 0, 

ot v qi [w ri + c) 


I dv u P2 ( A i 0 \ 

< — = a 2 Au - fj, 2 v + , C 4 - 1 - 2 ) 

ot V q2 W T2 


dw . u P3 

—— = a 3 Aw — UnW H-, x G SM > 0, 

^ Qt U® u ,r 3 


avec conditions aux bords de Neumann 


du dv dw 

dr) ’ drj drj 

et conditions initiales 


0 , 


x g r, t > o, 


u(z, 0 ) = ^(x) > 0 
u(x,0) = £ 2 (rr)>0 
w(x,0) = Ç 3 (x)>0 

où (j G C (fl) pour tout i — 1,2, 3. 

fl C est un domaine borné de classe C 1 de frontière <9 fl = T. 
Les constantes de diffusions: ai, 02 et «3 sont supposées positives, 
/ij, /i 2 et /i 3 sont des constantes positives, 

Pi, p 2 , P:i , çi, ç 2 , ? 3 , r'i, r 2 et r 3 : des indices non négatifs, 
u > 0 , c > 0 . 


(4.1.3) 


(4.1.4) 


4.1.1 Etude de l’existence globale en temps des solutions 


Présentons d’abord quelques hypothèses et conditions que nous allons utiliser dans la 
suite. 

Supposons que 


0 < pi — 1 < max 


<fp 2 min 


n \ 
\g 2 + 1’ r 2 ’ J 


, p 3 min 


/ r 1 qi 
\ r 3 + 1 ’ 



Posons 




I do , 

77= Pour z, j 

u a / d^dj 


1,2,3 


(4.1.5) 
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Soient a,(3 e t 7 des constantes positives telles que 

Pth + WhX 


a > 2 max < 1 , 


Pi 




a-1(3 + 1 a2 \ 1 lx - j. ,-r ± A 2 } > 


A 


12 


a — 1 7 + 1 




a — 1 


^23 — ^4l2^4l3 


(4.1.6) 

(4.1.7) 

(4.1.8) 


a (3 —j \ a 7 J \ ol 

Proposition 4.1.1 Le système(4-l-2) avec (4-3-3) et- (4-3-4) admet une solution classique 
locale unique ( u,v,w ) dans (0,T max ) x fl. Si T m ax < 00 alors 

lim (\\u ( t , OH^ + \\v (t, OH^ + ||w (t, .)IU = 00 

+/ J- max 

Preuve. Pour la preuve, voir Henry [18] et Ryan [45]. ■ 

Lemme 4.1.1 Soient p, q, r, s, m, et n des constantes telles que 


r > p — 1 
rm > (p — l)(n) 
rq > (p — l)(s + 1 ) 

Pour tous h, l, a, (3, 7 > 0, il existe C- = C(h,l,a, (3, 7 ) > 0 et 9 = 6 (a) G (0,1), telles 
que 


x- 


,p— l+CK 


a 


<(3~ 


X 


r+a 


yq+(3 £^+7 — r yS+l+(3 z n+^ 

Preuve. On a 

•p -1 


+ C 


T a \ 1 


ou 


x > 0,y > h, z > l. 


x* 


a- 


y q Z T 


= a 


x 


y s+l z n 


1 


(s + l)(p—1) _ (n)(p — l) 

y q y r z z r 


= 1(3 X 


y s+l z n 


P -1 


(s+l)(p-l) (n)(p-l) 

y q z r m 


= 1/3 X 


y s+1 z n 


P -1 


/ n X \ («+1)(P-1) _ (»)(P-1) _ 

B --- ü q Z r m 

y s+l z n ) 


pour un certain e qui satisfait 

( (s + l)(p-l) 

0 < e < min 


(n)(p — 1 ) 

+ q -LJLL - L + m 


\ 


(»+i) 


n 


,1 


p — 1 


\ 


(4.1.9) 
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Donc 


X‘ 


■p -1 


a -= a (BY 

y q z m 1 


P 


x 


y s+l z n 


P -1 


1_\ “ (y) (S+1) r (p - 1) -g+ 6 (s+l) ^Mfedi . m+eT1 



D’après (4.1.9) on trouve 


p — 1 


+ e < 1, 


(s + l)(p- 1) 


— Q + e(s + 1) < 0, 


et 


(n)(p- 1) 


— m + e(n) < 0, 


d’où 


(y) 


(s+DiP- 1 ) _ q+e{s+1) 


(s + l)(p-l) 

< h r 


— q+e(s+l) 

? 


(4.1.10) 


et 


P)(p-i) 

(Z) r 


— m+e(n) 


(n)(p—1) I / \ 
l ^ ---m+e(n) 


(4.1.11) 


D’autre part, nous avons 


et 




/3re 


/3re 

V > h => y 

a 

< h 

a 


7 re 


7 re 

z > l => z 

a 

< r 

a . 


(4.1.12) 


(4.1.13) 


(4.1.14) 
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Donc, d’après (4.1.10), (4.1.11), (4.1.12), (4.1.13) et (4.1.14) on aura 


x 


p -1 


a- 


y^z 1 


= a (13) r e \(3 


x 


yS+l z n 


r xv+e I]“ (y) ^-^ + l)^- m+En 



< a(/3) P r 6 ((3 


r \ ^++ e 


= a(PY 
< a (py 


(3 


(3 


y s+1 z n 


x 

y s+i z n 


yS+l z n 


J \ a h (s+iyp-1) _ q+<s+1) i („)(p-l) _ m+en 



X' 


p-1 


p-1 


+e 


+e 


1 \ “ (s+I-Kp 1) _ q+efs + 1 - ) (n)(p-l) . ÿre yre yre 

— | h L r y<xy<*z<*z°‘ 



x' 


1 \ “ ^ («+iKp-i) _ q+e(s+ i) t („)( P -i) 



•■/re _ (3re _"/re 


V p-1 , (a+l)(p-l) _ +e(a+1 )_gre ( n )(p-l) . 7 re / £ 

a (P)~~ h r a l LJ ^~ m+en -yr ( p 


I ' - m+en y* z «h«l- 

p-i 


y s+ 1 2 : n 


+e d y 13 z 1 




Posant 


p—1 ^ ( s +1)(p— il-Q-i-eCs-i-i)— ÊIH (n)(p— 1) . ryre 

C 1 — a (( 3 y^ h r ° i LJ y^-m+en-z-^ 


on obtient 


a- — < Ci ( y 

yq z m ~ 

Rappelons l’inégalité de Young-s 


p-i 


x 


y s+l z n 


p-1 


+ e /y0 z 7\ « 

x a 


fg < ef po + c £ g qo , avec /, g > 0 , £ > 0 , c e = e «- 1 ,p 0 , Qo > 0 , ^ ^ 


= 1 . 


Posons 

^ / = ( P 


cy 


X 


r \ yy+t 


ys+l z n 


9 = 


y 3 z 1 


x L 




On a 


- 1 -= 1 donne on = 

7)n <7n ^ u 


PO <?0 


1 - £=1 - e' 


Donc on obtient 

rP- 1 / 1 




Cl VV +1 W + U v 


_ p —1+re _re_ 

1 \ r —(p—1) —re / y/+ 7 \ / 1 _P^i_ £ 





alors 


X‘ 


.p-i 


a yi z m ~ v y s+1 z n 


x 


+ C 2 


x u 


y p zn) 


_re_ 

\ «( 1-^-0 


1 + 


OÙ C 2 = C*! 


p—1+re 
r—(p—1) —re 
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ce qui donne 


x- 


,p—l+a 


a- 


<p- 


X 


r+a 


+ c 2 


yq+Pzm +7 — ' yS+l+p% n +l 

Enfin, prenant e suffisamment petit on obtient 


( y 

W z V 


x- 


.p— 1+a 


a 


</? 


X 


r-\-a 


yq+Pz m + 7 — ^ yS+l+P 2: n +7 


+ c 


/ x a 


Lemme 4.1.2 Supposons que u (x ) > 0, v (x) > 0 et w (x) > 0. Alors pour tout ensemble 
d’indices p, q ,r, m, n, p et A satisfaisant X < p < m (non nécessairement positifs) et 
pour toute constante c> 0, on a 

—dx < c [ ^—dx + c -(p-a)\("*-p) [ ^—dx 

,q W r J n V n W p J n V ri VJ° 

où 

q (m — A) — n (p — A) r (m — A) — p (p — A) 

V = - i 6 = -• 

m — p m — p 

Preuve. Par application de l’inégalité de Young on achève la preuve. ■ 
Lemme 4.1.3 Soient u, v et w les solutions du problème (f.l.2)-(f.l.f). Alors on a: 



d_ 

dt 


1 ) 


2 ) 


vw 


< max 


< max 


1 

Ti 


o 


( 




1 


M2+M3 


92+^2+! ^2p 2 /lj2+»'2 + l 


— Ao. 


Preuve. 

1) On a pour k > 0 

— dx = —K 

U K 


u K+1 V dt 

1 


8 Aux 


= —K 


U K+1 \ 


( «i Au 


Pi u + 


u 


p 1 


v qi ( w ri + c) 


+ a ) dx 


= —nai (u K x ) (A u) dx + npi / — dx 


u 


■pi —«—1 


(4.1.15) 


(4.1.16) 


K 


V qi ( W ri + c) 


dx — KO 


u 


K-\-l 


dx. 
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Utilisant la formule de Green, on aura 


(u K x ) (A u) dx — (k + 1) [ ft . +9 \Vu\ dx , 


u 


alors 


1 [ —dx = 
dt / ■u'" 


-ac (ac + 1) ai 


1 ,2 . /' 1 


k+2 


U 


|Vd dx + AC/q / —dx 
/ u K 


u 


Pl-K-l 


— AC 


-dx — (J AC 


tt+1 


dx. 


u 


V 91 (w ri + c) 

Appliquant l’inégalité de Young sur le terme AC/p / Adx , on aura 


1 /"I 

ACyUi / — dx < ACCT / dx + AC<T -K p. , j i+1 1^1 ) 


CT 


U 


d’où 


- [ —dx < -K(K+l) ai 
dt / u K 


u 


P l-K-l 


m k+ 2 I V^l 2 dx + ACU J uK+1 

1 


dx + ACCT K /^ + |U| 


— AC 


U«1 (rU 1 + c) 


d x — a ac 


/■v+1 


dx 


u 


ce qui donne 


| 11* < w - >r i | fl |. 


Intégrant de 0 à, t . on trouve 


d 


Donc 


A / Adx ) dr < I acct K /r^ +1 |D| dr. 

CIL J u 


—dx < 2 max 
u K 


—dx , acct >* +1 |fi|i ) . 
si 


En élevant les deux membres à la puissance - et en faisant tendre ac —> +oo on aura 


1 

< max ( 

1 

— 

— 


1 

OO \ 

Ci 


CT 
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2) On a 


d 


dtj ( vw) 
n 


-(Kl = 


' —V (vw) 


fi 

= — ka 


vw 


s 2k 


-cin 


1 f 1 f U P2 

^ / yk+iy^jk d 2 I yky^jk I yq2-\-k-\-ly^T2-\-k 


(Kl 


— ka 


3 / ^ïîî A “’ + ^ 


v k w k 


-k 


w 


P 3 


'^; ( ?3+^'^; r 3+^+l 


<m. 


Appliquant la formule de Green pour les termes J vk + lyjk A vdx et J yk ^ k+1 A wdx ; 




A vdx = — J V 
n 


yk-\-\yjk 


Vvdx 


= ( fc + 1 ) f^K^k( Vv ) 2dx + k 


2jk-\-l^jjk-\-l 


VvVwdx. 


De la même façon on trouve 

f 1 


yjk-\-\yk 


Aiudx = (k + 1) 


1 9 

(Vw) dx + k 


n 

Donc 


yjk + 2yk 


yjk~\~\ yk~\~\ 


VvVwdx. 


d 

dt 


vw 


: dQ = — ka 2 (k + 1) 


-ka 3 (k + 1) 


yk-\-2yjk 


yjk~\~2 yk 


(Vu) 2 dx + k 


(Vw) 2 dx + k 


y k +1 yj k +1 


VvVwdx 


+k (/i 2 + n 3 ) 


w 


P2 


v k w k 


yQ2-\-k-\-l(jjjr2-\-k 


yjk-\-\ yk-\-\ 


dfl — k 


VvVwdx 


) 


w 


P 3 


yd3~fk yjT3-\-k-\-l 


(Ml 
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= -k(k + l)a 2 [^^{Vv) 2 dx-k 2 a 2 


yk-\- \ yjk-\-A. 


Vv'Vwdx 


~k (k + 1) a 3 / , (Vru) 2 cfe - A; 2 a 3 


yjk + 2yk 


yjk-\-\yk-\-\ 


Vv'Vwdx 


+k (/r 2 + /i 3 ) 


v k w k 


-k 


U 


P 2 


yQ.2 ~\~k -\-1 yj r 2 ~\~k 


d,Q — k 


u 


P 3 


^Ç3+/c^r3+A:+l 


rff2 


= / + J 


ou 


/ = —A; (A; + 1) < 2,2 


(Vu) cte - Ara 2 / ^ +lwfc+1 VuV W dx 


-A; (A; + 1) a 3 / (Vru) 2 cfe - k 2 a 3 


yj k +1 y k +1 


Vv'Vwdx, 


J — k (/x 2 + /^ 3 ) I k k 
n 

Pour /, on a 


u 


p 2 


V^w K J v q2+k+l w r 2 +k 

h n 


dkl — k 


U‘ 


P 3 


v 93+& w » - 3+A;+l 




/ = 


-A; (k + 1) 02 ^+ 2 ^ (Vu) - k 2 a 2 vk+1 1 wk+1 VvVw 
-k (k + 1) a 3 fc ^ a fc (Vru) 2 - A; 2 a 3 


^fc+l^yfc+1 
2 7 „ 2 „ 1 


Vu Vru 


— A; (k + 1) a 2 ru 2 (Vu) 2 — k 2 a 2 vwVvVw 


-k(k+l)a 3 v 2 (Vw) 2 -k 2 a 3 vwVvVw ) v ^ w 
dx 


k~\~2 q ! >k~\~2 


dx 


dx 


Q 


yk-\-2 yjk-\-2 


ou 


Q = k (k + 1) a 2 ru 2 (Vu) 2 + k 2 a 2 vwV (u) V (ru) 
+k (k + 1) a 3 v 2 (Vru) 2 + k 2 a 3 vwV (u) V (ru) 


= k (k + 1) a 2 [ru |Vu|] 2 + k 2 a 2 [ruV (u)] [uV (ru)] 
+A; 2 a 3 [ruV (u)] [uV (ru)] + k (k + 1) a 3 [u | Vru|] 2 
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Q est une forme quadratique par rapport à w | Vu| et v |Vw| , qui s’écrit sous la forme 
matricielle: 


Q = 


wV ( v ) 
vV ( w ) 



k {k + 1) Ü2 k 2 a,2 
k 2 a 3 k {k + 1) a 3 


wV [y) 
vV (w) 


Q est définie positive si, et seulement si, tous les déterminants principaux successifs Ai 
et A 2 de sa matrice des coefficients, sont positifs. 

Nous avons 

Ai = k (k + 1) a 2 > 0, 


A 2 


k (k + 1) a 2 k 2 a2 
k 2 as k (k + 1) a 3 
k 2 {2k + 1) a 2 a 3 > 0, 


donc Q est définie positive. 
Alors 


d 


dij (vw) 
n 


: dQ < J = k (p 2 + p 3 ) 


v k w k 


- k 


u 


P 2 


v q2+k+l w r 2 +k 


d,Q — k 


U‘ 


P 3 


v q3+k w r 3 +k +1 


dü. 


Utilisant le Lemme 4.1.2 pour le terme k (/i 2 + p 3 ) J ; 


nous avons 




-dx < c 


u 


u 


-dx + / -^—dx 

Jn v q w r J ü v n wP J Q v q w° 

d’après le principe de maximum il existe Cq qui dépend de IICiH^ , HC 2 II 00 de Il £3 II 00 
telle que u, w > C 0 > 0, d’où 
,p 


w 


-dx < c 
In v q w r J Q v n wP 


dx + c -(P-U\(--p) [ 


çyq+0 


On a ainsi 


1 1 

< 


u 


p 2 


V k W k /1 2 + p 3 J n V q2+k+1 W r2+k 


dx 


1 \ i2+ r 2+i r 1 

l-h + lh) Jn u 2 ^ 2 /^^ 1 


dx 
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avec 


P 

A 

V 

6 


0, q — k, r = k, c =-, m = p 2 , n = g 2 + k + 1, p = 72 + k 

S2 + P3 


ç 2 + r 2 + 1 

q (m - A) - n (p - A) .92 + 1 

-= k + A- 

m — p p 2 

r (m - A) - p (p - A) _ ^ ^ 

m — p P 2 


alors 

2k -, 

d f 1 / 1 \ «2+' - 2+ 1 r l 

dtj (vw) kdn ~ * U 2 + / 4 J Jn u 2kp i/v + n +l dx 

n 


d’où 


d 

dt 


dfl < k 


(vw 




^2 + p3 


Q2+ r 2 + 1 


-1 


|fi| 


2fcf>2 / Q2 +r-2+1 


Intégrant de 0 à, t. élevant à la puissance l et faisant tendre k vers + 00 , on trouve 


1 

vw 


< nrax 

OO 


^2^3 


l l 2 + p3 


12 + r 2 + 1 


Aî 


2p2/q2+r2+l 


A-2- 


Lemme 4.1.4 (voir [36]) 

Soit L + (0, T) l’ensemble de toutes fonctions intégrables f (t) > 0 dans (0 ,T) telles que 
la quantité 

K[f}= sup [ e~A T - T lf (r)drr 
0<t<T J 0 


est, finie. 

Supposons 0 < 0j < 1 et 7 ^ G L + (0 ,T) (j = 1J). Soit W = W (t) une fonction 
positive dans [0,T) satisfaisant l’inégalité différentielle: 


^ < -fiW ( t ) + Y jlj (t) W°> (t), 0 < t < T. 

3 =1 


Alors 

W (t) < k, 0 < t < T, 
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où k est la racine positive de l’équation algébrique en x: 


x é K b il x0j - W (°) 

3 =1 

Théorème 4.1.1 Soit la fonctionnelle 

f u a 

Lit) = / — a — dx. 

y ’ J vdw'i 

ü 

associée au système (4-1-2). 

Sous la condition (4-1-5) , la fonctionnelle Lift ) est uniformément bornée dans l’intervalle 
[0,T*], T* < T max , où T max (IKIL , IK ; o||oo , Hwolloo) définit le temps éventuel d’explosion. 


Corollaire 4.1.1 Sous les hypothèses du Théorème 4-1-1, la solution du problème (4-1-2) 
est globale et uniformément bornée dans Ox ( 0 , oo) pour des conditions initiales positives 
dans C (fl) . 


Preuve de théorème 4.1.1. 

Dérivant L(t) par rapport à t, on obtient 


d 


1-Mt) = 37 


d ( 


dt 


dt 


i—l 

fvhw't J 


dit 


v 2/3 w 2 'r 

u a ~ l ( du 


- fiu^^w 1 ^- - 7 ) dLl 
dt y dt 1 ^ 


a 


vhw'i V dt 


3i7 ~fi 


U 


f dv\ u° 

I eu. J h . 


vh +1 w~i \dt J v> 6 un +1 \dt J 


dw\ 


dt 

) I dtt 


( u a ~ l 


a 


vhw^ \ 


( a± Au 


Pi u + 


U 


P 1 


v ( n (w r 1 + c) 


V 


-fi 

-fi 

-7 


u 


v / 3+! w 7 

U a 

yh +1 w 7 
u a 


a 2 Av — p 2 v + 

> 

a -2 Av — p 2 v + 


u 


P 2 


+ a 

v q2 w r 2 ) 


\ 


u 


P2 


V q2 W r2 7 


V /9 W 7+1 V 


Uaw - + ff-A 


dLl 
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donc 


d 

dt 


Lit) 


( u a 1 u 

a±a— — : Am — /i x a 


yPw' 1 


yPw' 1 


- 1■ Ol~ 


u 


,pi+a—l 


+ a a 


u 


a —1 


v<n+Pw' y (w ri + C / 
a2 ^ v@ +1 w' r ^ V ^ v q2+ @ +1 w r2 +')' 

“^^7+!^ + ^ y Pyj'f ~ 7 n®+/3 u ,r 3 +7+l 
a-1 

«2/3 ^ „ Am - a 3 7 


yPw'i 


dü 


u 


ai a—s—A m 
V /3 W 7 


—/i x a 


M 


v / 3+! w 7 

„ m l 
+ t l 2P~ 


H-a 


u 


Pl+a—l 


\ ÎJ« 1 +%'l'(w ri + c) 


V /3 W 7 ' r2r V P W l 

v P2+a 

-P- 


+ /^s7 


vPw"1 +l 

u 

V&W 1 


Aiu 


v q2+0+l w r 2 +'y 

= /+ J 


-7 


u 


P3+Ot 


yq3+0yji'3+'y+l 


+ aa 


U 


a —1 


dn 


vPwl J 


L -X »_/ l^C/U U 

Y 


a 


J 


fu a ~ l f u a 

aia ] AudSl ~ 02/5 J 

n n 

f u a f u a f u a 

- ,ha J in + ^J ^ da + lhl J dS1 

n 

/ u P2+a f u P3+a 

v q2+P+l w r 2 +'Y^ 7 J v g 3 +/3 u ,r 3 +7+l 

n n 


n 


f u tn+a ~ l 

a J U «i+/ 9 l y 7 ( lü ri _|_ c) 

n 

Alors 

r^cx—l r yCx r yCx 

^ Audü - a 2 /J J i - « 37 J ^^AwdSl, 

n n 


-1 


M P3+« f u a-l 

5 - ——dit + aa / —5 - 

3 w r-3+7+l J yPyjl 

n 


u a ^ 

dû. 


aia / —7—/ 

n 


n 


J = (-/^a + p 2 /3 + // 3 7) L(t ) 


-|-a 


u 


pi+a-l 


v < n+Pw' 1 {w ri + c) 


dfl — /3 


u 


P2+Ot 


X;92+/3+l u ,r 2 +7 


dO — 7 


U‘ 


P3+OI 


yq3+P w r 3 +^+l 


d,Q + aa 


M 


a— 1 


y0W^ 


dit. 


En appliquant la formule de Green, tenant compte des conditions aux bords de Neu¬ 
mann, on aura 
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v a-l /' / 7 .a-l 

-Awdfi = - / V f ) V («) dfi 


V^W 1 



v^w 1 


n 


-/v(«: 


(a — 1) m q 2 v^w 1 'V ( u ) 
v 2 i 3 w 2 ^ l — (j3vP~ 1 w' y 'V (v) + 7 v /3 w 7_1 V (iü)) 


/ (“- 1 ^i v “i 2 - / 3 ^ v w v w X 

m"- 1 


V 


7 fl 7+1 V (u) V (w) 


dVL 


dn 


donc 


u 


a—1 


v^w 1 


Audi} = — 


«“- 1 


V 


-7 


n /3 w 7+1 


V (u) V (w) 


dfi. 


^A««=-yv(,)v(^) 


dO 


■-- - / V (w) 


i «,i / (/? + 1) v^w^V (v) 

2/9+2 2'y i (u) - ( I u a d,Q 

y2p-\-2yj2^( I 


/ ^ 7/ Q! 

‘ V(«)V(n)-(/3 + l)^^|Vn| 


+7n /3+1 w 7 (w) 

2 \ 



vfi+ 1 W 1 


~1 




•y/3+ 1 'U;7+ 1 


V (v) V (w) 


dfi. 


1 


Ainsi 


u 


v /3+ 1 U /7 


Audi} = — 




|v»|- 


-i 


U 


vP +1 w7 +1 


V (v) V (w) 


dQ. 


n 


* 7 i a 

-^—Awdtt 

v /3 u ,7+ 1 


n 


v <“ )v (+Sî) 


An 
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vh 


au a l v^w 1+l \/ (u) 


v 2/3 w 2 7 +2 1 ^ W 7+1 + V H+ 1 ) V P) 


dn 


U 


( u a 1 


a 


vPw' 1+l 




u a \ 

v(“) v W-^lîvïïVWvW ' 

- (T + !) |V '" 12 


dn. 


Par suite, 


/ yCV 1 ^.O! \ 

I T—r / \ T- 7 / \ n T- 7 / \ T- 7 / \ \ 


U 


V0w' y+1 


A wdn = 


a 


V^W"! +l 




v(“)vW-^ +lro1+1 

~ < 7 + 1) Hï |Vm|2 


v (v) v H 


dn. 


Alors on obtient 



-a\a I | (a 

h 

o u °~ l 

+&2P 



d i v “i 2 - 7 (M) v (t,) - (m) v w ) dÇl 

v («) V (v) -1/3 + 1) Hï l Vi ’| 2 - H v (t,) v W ) da 
n 

p / (x — 1 (X OL \ 

+<i37 J v (u) v (m) _ /? H^îî v w v w - h+1) 1 v " ,|2 J dî! 

n 

/ -aja (« - 1) 'L— | Vu| 2 + a/3 (a, + a 2 ) V (n) V (n) \ 

r U J GL ~ ^ 

+« 7 (a, + a 3 ) V (■ u ) V H 

^ W + 1) ^1 Vv l 2 - ^ ^ + fl 3) V ^ V (w) 

U a 2 

\ -W(7 + 1) ÏSS ^|VH y 


O 


dû 


f — a\a ( a — 1) v 2 w 2 | Vw | 2 + a/3 (ai + a 2 ) vuw 2 V (u ) V (v) ^ 
+«7 (ai + a 3 ) v 2 uwV ( u ) V (w) — a 2 /3 (f3 + 1) n 2 w 2 | Vn | 2 
— 7 P (a 2 + a 3 ) vu 2 wV (v) V (w) 

V -a 3 7 (7 + 1) H 2 |Vw | 2 / 


u 


a—2 


yd+2yj'Y+2 


dn 
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( — aia (a — 1) [vw |Vu |] 2 + a/3 (ai + 02 ) [vw'V (u)] [uwV (r>)] ^ 
+«7 (ai + a 3 ) [vw'V («)] [uvV (ia)] — a 2 /3 (/3 + 1) [uw |Vi :|] 2 
-7 (3 (a 2 + a 3 ) [uwV (u)] [uvV (w)] 

\ -«37 (7 + 1) [uv | Vw |] 2 / 


u' 


a .-2 


i;/3+2^y7+2 


dfi, 


= ~ Q 


u' 


a-2 


i;/3+2yj7+2 


dfi, 


ou 


^ aia (a — 1) [îw | Vu |] 2 — a/3 (ai + a 2 ) [iwV (u)] [muV (u)] ^ 
—«7 (ai + a 3 ) [vwV (u)] [uvV (w)] + a 2 /3 (f3 + 1) [uw |Vi ;|] 2 
+ 7 (3 (a 2 + a 3 ) [uwV (u)] [uvV (w)] 

\ +a 3 7 (7 + 1) [uv | Vw |] 2 / 


Q est une forme quadratique par rapport à: vw \Vu\ , uw |Vu| et uv |Vw| , qui s’écrit 
sous la forme matricielle: 



^vwV (u)^ 

T 

( 

f ai a(a 1) a(3 ai+ 2 a2 a 7 ai + a3 ^ 


^vwV (u)^ 

\ 

Q = 

uwV (v) 



a(3 ai + a2 a 2 /3 ((3 + 1) /3 7 a2 + a3 


uwV ( v ) 



yUvV ( w ) j 


V 

^ -« 7 ^ /?7 ^ «37 (7 + 1 )) 


yUvV ( w ) j 

7 


Q est définie positive si, et seulement si, tous les déterminants principaux successifs 
Ai, A 2 et A 3 de sa mutrice des coefficients, sont positifs, 
où 


Ai = ai» (a — 1 ), 


Ai > 0 d’après (4.1.6) . 


Ao — 


101+02 


aia (a — 1) — a/3 2 

a 2 (3 {(3 + 1) 
(a — 1) {(3 + 1) 


= a 2 (3 2 aia 2 


a 


fi 


- A 2 

Ji l2 
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d’après (4.1.7) nous avons A 2 > 0 . 


A;» — 


oio (o — 1) — af 3 ai ~^ a2 
-o/ü^+jp 1 a 2 (3 (f3 + 1) 


- 07 ^^ 




h q ^r A - 

037 (7 + 1) 


aict (o — 1) 


a 3 o (0 — 1) —o/ 3 ai )p 2 


OiO (0 — 1) — 07 

o/ 3 ai + a2 a 2) d (/3 + 1) 

1 ^ . 1 

1 

o/ 3 ai +“ 2 /57 a2 + as 

t 


a\a (o — 1) —cr/3 a - 1 + ct2 oio (o — 1) 


- 07 ^^ 


/?7 


«2+a3 




-07^+^ 037(7 + 1) 


o 


2 n2 2 

0/57 010203 


(a - 1) 


ot — 1 / 3+1 y^2 j | a—1 7~l~l ^2 


Q; 7 


a /3 

_ (+T^23 — A 12 A 13 ) 


L 13 


et par (4.1.8) et (4.1.6), on obtient que A 3 > 0. 

Estimation de J 


J = (-/+o + p 2 /3 + p 3 7 ) L(t) 


+0 


u 


P\+OL—l 


Vl 1+ Pw’ y (w ri + c) 




U 


P2+a 


v q2+P+l w r2+'y 


dû — 7 


ÎJ- 


P3+« 


1)53+/% r 3 +7+1 


df2 + (Ta 




a—1 




dû 


n n n fi 

Le principe de maximum nous fournit l’éxistence de C 0 qui dépend de H^iH^ , ||£ 2 | 
et de H^IL telle que v,w > C 0 > 0, donc on a 

o — l (3 7 o — l /3 + 7 

u a ~ x ( U a \~ fl\a fl\â / ( / 1 N 

VPW'T ~ J [vJ [wJ ~ 

ou encore 

o — l /3 + 7 

ou 


1 \ O 

C 0 J 


<c,/A!_Y 


v^w 1 


\vPw' y / 


o 


C-3 = ( VV 
c-n 


o 
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On a 


J = {-^a + h 2 /3 + ntf) L(t ) 


+a 


~1 


u 


pi+a— 1 


'y<îi+/3'U;7(-u; ri -f- c) 

,P3+a 


dQ — (3 


U 


P2+a 


v q2+P+l w r2+l 


dn 


w 


-y93+/3 w r3+7+l 


df 2 + (xa 


U 


a—1 


v^w 1 


dn. 


n n 

En appliquant le Lemme 4.1.1 on obtient 


pi+a—l 


pi+a—1 


P2+a 


a 


f U° Y 

v qi+ dw' ï (w ri + c) — a v qi +Pwi +ri ~ v q 2 + P + 1 w r2+ ' y + C \vPwi J 


donc 

J < (-/x 1 a + n 2 p + /x 3 7) L(t) + /3 


u 


P2+a 


v q2+p+l w r 2 +l 


dn+ c 


{Y)’ da 


a—l 


-\-(X<J / C 3 


f—Y 

\ vd'Uf’ J 


a 


dn — /3 


u 


P2+a 


yq2+P+lyjr2+'f 


dn — 7 


U 


P3+« 


^<73+/3 w r 3 +7+ 1 


dn, 


il s’ensuit 


J < (—+ h 2 (3 + P37) L(t) + C 


YY) da+a Y C3 


(—) 

yvdw'y J 


a — l 


a 


dn. 


Compte tenu de l’inégalité de Hôlder, on obtient 

r / \ e 

f u a y 


1-0 


C 


K**) m ^\iYY) ia 


(c)^ dn 


donc 


f / a \ 8 

c ( ïYY ) dn ~ C ±YY )) 6 , où <74 = 0101 


1-0 


On a 


a 


(J—') 

\vPw7 J 


a — l 


a — l 


a 


dn < 


(—) 

\ydw^ J 


a 


d.n 


(c 3 ) a dn 
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alors 

a a ~ 1 a - 1 ^ 

j C3 {j^J a dn < C *W))~ ’ où C 5 = C 3 |fl|« . 

n 

On obtient ainsi 


a — 1 

J < (—/^a + [ 12/3 + p. 37 ) L{t) + C 4 (A(t)) + otuC^ ( L(t )) 


ce qui donne 


a — 1 ' 


J < (— fi-yCx + p 2 (3 + P 37 ) L{t ) + Cg [ (ü(t)) + au (A(t)) 


a 


Sous les conditions (4.1.6), (4.1.7) et (4.1.8) on aura: 


a — 1 ' 


7/ (t) < (—/^a + /i 2 [3 + /r 3 7) L(t) + ( 7*6 ( L(t )) + crcr ( L(t )) 


a 


et comme — /x x a + fj, 2 P + /i 3 7 < 0 d’après (4.1.6), il en résulte 

L'(t) < C 6 L 0 (t) + CgacrL^(t). 


(4.1.17) 


En appliquant le Lemme 4.1.4 sur (4.1.17) on trouve que L{t ) est bornée dans (0, T max ); 
i.e. L(i) < 73 , où 73 dépend de la norme-L 00 de £ 1; £ 2 et de £ 3 . ■ 

Preuve du Corollaire 4.1.1. 

Supposons Au = dAu — /qu, alors A est le générateur infinitésimal d’un semi-groupe 
holomorphe {e tA |0 < t < 00 } dans L a (fl), de domaine 

D(A) = iu \u G W 2 ’ a (fl) , ^ = 0 j , 

où a = N (N + 1) /2 et N est la dimension de R N . Par ailleurs, il existe des constantes 
ç > 1 et 5 > 0, telles que 


Ah tA 


< ç ( 1 + t 


St 


t > 0 , 


(4.1.18) 


où ||.|| dénote la norme de L a (fl). Ainsi la première équation du système (4.1.2) avec 
du 

— — 0 et u (x, 0 ) = ^i(a;), est équivalente à l’équation intégrale 


u{t) 



u pi (s) 

v qi (s) ( w ri (s) + c) 


+ a 


ds. 


0 



Utilisant le Théorème 1.6.1 cité dans [18] et (4.1.18), on obtient 

nia xu(.,t) < A^u{t) 
n 


(4.1.19) 


+ fç [l + (t-s)-^e- s ^ s \ 


u pi (s) 


v qi (s) w ri (s) 


+ O 


ds. 


Supposant m > ap et rriq > pn et rm > kp et appliquant le Lemme 4.1.2, on trouve 


u 


pi 


v qi w ri 


U 1 


pi a 


v c n a w ria 


(Kl 


n 

< c 

< / 


U 


(Pia-A) r u 


A 


<Q v n w k 


dx + c (™-pi“) 


v q w e 


dx 


u 


l n v n w k 


dx + 


v q w e 


dx 


ou 


aqim — anp\ ar^m — akp\ 

p = - , 0 = -, A = 0. 

m — api m — ap\ 

Grâce au Théorème 4.1.1 il existe q 4 qui dépend de la norme-L°° de £ 1; £ 2 e * de £ 3 ,telle 


que 

D’après (4.1.16) on a 
/• 1 


u 


! n v n w k 


dx < 7 4 . 


v q w e 


dx < 


vw 


Alors 


donc 


u 


pi 


v qi w ri 


< A 2 tq t) est constante. 


< (73 + A t) “ = : A 4 


maxu (., t) < ç ^1 + t 2 j e 5 t ||£i||+çe St f e Ss .A^ds + çe 5 t fe ôs .ads (4.1.20) 


Le Théorème 4.1.1, inégalité de Hôlder, (4.1.16) et (4.1.20) impliquent 

max-u (., t) < M ^1 + t . (4.1.21) 

Il résulte de (4.1.21) et du principe du maximum que la solution est globale et unifor¬ 
mément bornée dans Ox (0, 00 ). ■ 
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4.1.2 Explosion en temps fini des solutions pour le système 


EDO associé 


S'il existe un temps t max < oo, tel que les solutions sont bien définies pour tout 0 < 
t ^ Tmax, ou 

lim ||u (t ,.) + v (t , .) + w (t, 0 IL = oo, 

J- max 

alors on dira que les solutions explosent en temps fini. 

Supposons dans notre problème (4.1.2)-(4.1.4) que c = 0. 

Considérons d’abord le lemme suivant déduit du principe du maximum (voir [43] 
page 297): 


Lemme 4.1.5 Soit (u(t,.),v(t,-),w(t,-)) la solution de (4-1.2)-(4-1.4)■ Alors pour tout 
( t , x ) dans ( 0 , T max ) x fl, on a 


u ( t , x) 

> 

e _Ml< min (^i(x)) > 0 , 

v ( t , x) 

> 

e - ^ 2 * min (£ 2 (x)) > 0 , 

w ( t , x) 

> 

e - ^ 3 * min (£ 3 (x)) > 0 . 


(4.1.22) 


Théorème 4.1.2 Supposons 

pi -1 min r_ (pri _ 1 

max (p 2 ,P3) L min (? 2 , 93 , r 2 , r 3 ) - max ( 91 , ry) + gr/y ’ J ’ 

alors il existe un espace-indépendant des conditions initiales uo,vq,wq, tel que la solution 
(■ u,v,w ) = (u (t ), v (t ), w (t)) du problème (4-1.2)-(4-1.4) explose en temps fini. 


Preuve. Considérons l’espace-indépendant des solutions de (4.1.2)-(4.1.4), i.e. les 
solutions de système EDO associé sans diffusion. Pour des conditions homogènes spa¬ 
tiales u 0 ,v 0 ,w 0 > 1, supposons pour contradiction que T max (u 0 ,vo,w 0 ) > 1- Dans ce 
qui suit, toutes les constantes positives C 7 , Cs,... sont indépendantes de uo,vq,wq. 
Pour a, /3 ,7 fixés non négatifs, soit A = ap 1 — f3p 2 ~ IP 3 et L (t) = 

Nous avons trouvés précédemment que 

/ A o A “"a \ 

I n, rv -A 11 — n~ H -/\7; — ^ /\ nn \ 




m = 


u u 

aia —;—Au — aoS—T —— 

u „ u' 


«37 


—/i x a 


Soi 


u 


pi+a—l 


P 


V^W 1 
U 


S p 2 P 

P2+OL 


yQ.l+Pyjl+'f'l yQ2+P+^yjr2-\-J ^Ç3+/3^r3+7+l 


vPw 1 
7 


+ /D7 


ff% 7+1 

u 


A w 


u 


P3+CK 


+ aa 


u 


a—1 
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/ 


U 


,pi+a—l 


\ yQi+Pyyl+ri 


-fliCX 

~P 


U 

V^W 1 

U P2+a 


+ V>2p 


U 


U 


V&W 1 
U P3+a 


\ 


U 


a —1 


donc 


ïT7 + ft > 7 

U/ Uj U/ I 

- 7 -h (T CX - / 

y12+P+Xyjr2+l yQ3+/3yjr3+-i+l yP'W'Y / 


/ 


U 


U 


L'(t) + A L(t) = 


~HiCt 

■P 


, u^w 7 „ w v w 

M Pi+o-l yP2+a yP3+a U a ~ l 

+0;- /3 -'Y-f- cra;- 

V 9i+/3 w 7+n t) Q2+/3+l u ,r-2+7 ' v <?3+/3 w P3+7+1 

Xi 0 Zi a 

+°Al -â-Z - Pu2~â~Z - 7/^3 -S-Z 




U 


yPw^ 

u P3+a 


\ 


L A % 7 


U 


.pi-l+a 


= a 


? ,P2+a 

-Æ-7 

'))Çi+/3'j/)7'i+7 ^'?i92+/3+1,/,)’2+7 ' 


V&W 1 

u P3+a 


V^W 1 

u a ~ l 


\ 


) 


Envisageons deux cas séparés: 

Cas 1 : pi — 1 > rnax (p2, P:s) ■ Appliquant (4.1.23) avec a = 1, et prenant /3 ,7 suff¬ 
isamment grands , utilisant (4.1.22) , £ 2 > A 3 A 1 et l’inégalité de Young; nous obtenons 
pour tout t G [0,1] 


a u Pl 
3 v'?i +/% î7 + 7 


a ( u P2+1 

3 \^yQ2+l+/3yjr2+l+'ï 


P 1 /P 2 +I 


v k w h > P 


U P2 + 1 

v q2+P+i w r2+'y 


-Ch 


où k = (q 2 + 1 + P) ^7 - qi - P > 0 , h = (r 2 + 1 + 7 ) ^ - ri - 7 > 0 

. . 1_El_ pi 

(J 7 = (Aj P1-(P2+!) ^Pi-(P2+!)_ 

Avec les mêmes arguments nous obtenons 

a u Pl ^ u P3+1 

3 y<n+Pw ri+ ' y ~~ ^W r 3 + 7 + 1 r;'?3+^ 8 

, -1_pi_ pi 

OÙ Cÿ = (f) Pi-^s+i) 'ypi —(p 3 +l) . 

D’autre part, en utilisant (4.1.22), on arrive à 

a u Pl _ a / u \pi > / u \pi 

3 pî i +%n+ 7 3 ViAw 7 / ~ 9 VrAw 7 / 


où m = (pi — 1) P — qi > 0 , n — {pi — 1) 7 — 17 > 0. Donc on a 

U pi ^ u P2+<* ~ U P3+* / U \Pi 

Q y<îi+^ w Pl+7 — ■y<î2+/3+l U ;P2+7 ^ ^ w r 3+l+ 1 y q 3+0 9 Xypyj'ï ) 

de (4.1.23), on obtient 


L'(t) + A L(t) > C 9 ( L(t)) Pl -C 7 - C 8 , 
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ou encore 


L\t) + A L(t) > C 9 ( L{t)) pi - C, 

où C = C 7 + C 8 . 

Prenant L(0) suffisamment grand, ce qui implique l’explosion de a avant t = 1, ce qui 
contredit avec T rnax > 1 

Cas 2 : 

Pi -1 < max (p 2 ,p 3 ) et (jh - 1 ) [min (q 2 , q 3 ,r 2 , r 3 ) - max (q t , ri) + gpri] > q 1 r 1 max (p 2 , p 3 ) . 
Affirmons qu’il existe des constantes Cio, Cu > 0 telles que, si 

*max(p 2 iP3 )~Pi +1 v. r t m in(g2,<?3)-9l t min(r 2 ,r 3 )-n 
SI A C 10S2 S3 » 

alors 

M max(p 2 ,P3)-Pi + l > CiiU m in( 9 2>5 , 3)-9l u , min ( r 2,r 3 )-ri 0 < t < 1 


(4.1.24) 


(4.1.25) 


Pour prouver cette affirmation, prenons Z ( t ) = e M L ( t ) et appliquons (4.1.23) avec 
a = max (p 2 ,p 3 ) - pi + 1 > 0 , (3 = min (q 2 , q 3 ) - q 1 > 0 et 7 = min (r 2 , r 3 ) - r\ > 0 , 
alors, pour tout t G [ 0 , 1 ] , 

Z' (t) < 0 ==► Il ( t) + XL ( t) < 0, 


ce qui entraine 


M max(p 2 ,P 3 )+a 

v mm(q2,q3)+f3 w tnin(r2,r 3 )+'y 


> 


u Pi-l+a 


d’où 

Z(t) > > e-' A 'c(/? + 7)“' =: Cio- 

Par conséquent, on a Z (t) > Cio dans [0,1], donc l’affirmation est prouvée en prenant 
Cu = e^Cio. 

Maintenant, supposons (4.1.24). Utilisant la première équation dans (4.1.2) et (4.1.25), 
nous déduisons 


u 1 + pi u > 


u 


pi 


v qi w ri 


DI 31 Tl max(p 2 ,P 3 )-Pl + 1 „ 

> Cl 2 « y 1 min( î 2.î3' r 2. r 3)-max( ïl ,r 1 ) — C\ 2 U°, 


0 <t<l 
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où 

_ nn. _ 

/~1 _ / ^fmin(q 2 ,q , 3,’’2,'r3)-max((j 1 ,r 1 ) 

12 — O 11 

max(p 2 ,p 3 ) ~Pi + 1 

V = Pi-qiri—, -\- 7 - 7 

mm (q 2 , q 3 , r 2 , r 3 ) - max (<?i, ri) 

1 (pi - 1 ) [min (g 2 , ? 3 , r 2 , r 3 ) - max (gi, n) + gyri] - gin max (p 2 , p 3 ) > 1 

min (g 2 , g 3 , r 2 , r 3 ) - max (gi, n ) 

En prenant suffisamment grand, on obtient l’explosion de a avant t — 1, ce qui 
contredit avec T rnax >1. ■ 

Remarque 4.1.1 Le modèle biologique de Phyllotaxie (4-1-1) est un cas spécial du sys¬ 
tème (4-1-2) pour 

p 1 =2,q 1 = l,r 1 = l,p 2 = 2,q 2 =0 ,r 2 = 0 ,p 3 = l,ç 3 = 0 ,r 3 = 0 . 

Il est clair que ces constantes réalisent les conditions (4-1-5) de l’existence globale. 

4.2 Existence globale des solutions d’un système de 
type Gierer-Meinhardt couplé 

Selon le modèle d’origine de Gierer-Meinhardt à deux équations (3.3.2), Crystal Cooper 
[ 8 ] a proposé un autre modèle théorique composé de quatre équations de Réaction- 
Diffusion de type Gierer-Meinhardt couplées avec une modification de Conway-Cooper, 
où on ajoute le lien de la différence entre la première et la deuxième équation à la 
troisième comme suit: 

d [rti] /dt = d Ul d 2 [ui] /dx 2 - n [tp] + c 0 p [uif / [u 2 ] + p 0 p, 
d [u 2 \ /dt = d U2 d 2 [u 2 \ /dx 2 - v [u 2 \ + dp' [u ^ 2 , 

d [u 3 \ /dt = d U3 d 2 [u 3 ] /dx 2 - p [u 3 \ + c 0 p [u 3 ] 2 / [u 4 \ (4.2.26) 

+PoP + ky([ui] - [u 2 \ + floor ), 
d [u 4 \ /dt = d U4 d 2 [u 4 \ /dx 2 - v [u 4 ] + dp' [u 3 ] 2 , 

où 

[mi] , [it 2 ] i [ u :s] et [?i 4 ] représentent les concentrations de quatre substances u 4 , u 2 , u 3 et 
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U 4 . 


p et p' sont les concentrations de source de a\. u 2 et u 3 . « 4 , respectivement. 

Le terme p 0 est une constante qui représente la production basique de a \ et u 3 . 

Le symbole p représente la destruction ou le déplacement de u 4 et u 3 , et v représente 
la destruction de w 2 et u 4 . 

co et d sont des constantes de production de a\. u 2 et u 3 , u 4 . respectivement. 
d ui , d U2 , d U3 et d Ui sont les constantes de diffusion de v,i , u 2 , u 3 et u 4 , respectivement. 
ky([ui] — [ 112 ] + floor) représente le lien de la différence entre la première et la deuxième 
équation où ky et floor sont des constantes contrôlent la taille de cette différence. 

Considérons maintenant le système de Réaction-Diffusion de type Gierer-Meinhardt 
qui généralise le modèle (4.2.26) pour étudier l'existence globale des solutions: 


,p 1 


( dui . , u\ 

= aiAwi - pui + b 1 —z- + a, 


dt 


du 2 

~dt 


u 


a r U 1 

= a 2 Au 2 - vu 2 + 02 ^ 7 , 


r 1 


U 


,P2 


d U 3 Un 

— = a 3 /\u 3 - pu 3 + b 3 —^ + cr + a(ui - u 2 + b), 

Ui U ^ 


ÔU 4 Ul 2 

— = a 4 Au 4 - uu 4 + 64 — , 

u~C 

avec conditions aux bords de Neumann 

dut d u 2 du 3 d u 4 

= 0 , —- = 0 , = 0 et 


x G O, t > 0, 


x g n, t > 0, 


x G O, t > 0 , 


x g n, t > 0 , 


(4.2.27) 


dp dp 

et conditions initiales 


dp 


dp 


= 0, 


x g r, t > 0 , 


(4.2.28) 


ui (x, 0 ) 
u 2 (x, 0 ) 
u 3 {x, 0 ) 
u 4 ( x , 0 ) 


<Pi( x ) > 0 , 
(p 2 (x) > 0 , 
(p 3 (x) > 0 , 
(p 4 (x) > 0 , 


ieO, 
x e O, 
ieO, 
igO, 


(4.2.29) 


ou 


fl C est un domaine borné de frontière T régulière. 
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Les constantes de diffusions: 01 , 02,03 et o 4 sont supposées positives. 

/i,z/ et a sont des constantes positives, 

Pi, p 2 , q\ , q 2 , ri, r 2 , Si et s 2 : des indices non négatifs avec pi > 1, p 2 > 1- 
61, 6 2 , 63, et 64 des constantes positives telles que b\ < b 2 et 63 < 64. 

4.2.1 Etude de l’existence globale 

On va démontrer que si: 


ri > Pi - 1, r 2 > p 2 - 1, riçi > (pi - l)(si + 1) et r 2 ç 2 > (p2 - l)(s 2 + 1), (4.2.30) 


à l’aide des résultats du chapitre II, le problème (4.2.27)-(4.2.29) admet des solutions 
globales bornées en tout le temps. 

Lemme 4.2.1 (Idem, que Lemme 3.3.1) Soit p,q,r et s des constantes telles que, r > 
p — 1 et rq > (p — l)(s + 1). Pour tous h, a, (3 > 0, il existe c = c(h,a,/3 ) > 0 et 
9 = 6 (a) G (0,1), telles que, 


u 


p—l+a 


a- 


V' 




<p 


V 


u r+a fu 
+ C ( —ô 


)^+l+/3 


V 


ou 


u > 0 , v > h. 


Théorème 4.2.1 Soit la fonctionnelle 


U (i)= l n fd X+ f'fdx. 


n 


u 


n 


Ua 


associée au système (4-2.27) avec 


1 Pi v> \ f (3 2 v 

a 1 > 2 max 1 ,- f ,a 2 > 2 max < 1 , 


p 


/i 


1 ^ / Oi + 02 \ , 1 

~ô~ ^ l ô / _ _ ) _ b ~ô 


Pi \ 2 ^/âfâ) J 


( 03 + 04 V 
P2 " v 2^/03O4 J 


— > 


alors pour 


Pi - 1 


< min 




1 \ P 2 1 . ( q-2 - 

,11,- < mm | -r, 1 


ri \si + l’ J’ r 2 . \s 2 + l : 

la fonctionnelle u(t) est uniformément bornée dans l’intervalle [0,T*], T* < T n 

où T max ( || (Pi || 00 , Il <P 2 II 00 ji II P 3 II 00 > Il Pi II 00 ) définit, le temps éventuel d’explosion. 


(4.2.31) 

(4.2.32) 

(4.2.33) 
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Corollaire 4.2.1 Sous les hypothèses du Théorème 4-2.1, la solution du problème (4-2.27)- 
(4-2.29) est globale et uniformément bornée dans Ox (0, oo) pour des conditions initiales 
positives dans C (fl) . 


Preuve de théorème 4.2.1. 


üü(t) — I ^2—dfl+ I ^-diï = Ui(t) + üü 2 (t), 

1 V 


OC1 


U 


ai 


fl 


fl 


P 2 


OU 


f u ai r u a2 

u)i(t) = I -j^dVL et cu 2 (t) = / -|^o?fl. 


u 


U ,, 


Dérivant u(t) par rapport à t , on obtient 


ou 


jt- Il (if J dn+ Il (ü • dn - 1 + 

n y / fi y 


7 = 

n x 2 


J -/s ï 1 "- 

n v 4 




n 


dt [ uî 1 


clfl 


/^ (“‘““"A (^) ^ (^)) <*«• 

n 2 


u 


ai —1 


«1 - 


U 


A 


( 

aiAui 


u pi N 

■ +A + 61-57 + 0- 
«9 


/ 

^1 /3j + l ( 

u 2 V 


VU 2 + 62^7 dD 


,n ’ 


u 


( 


fl 


< 1_1 A O < A /O t < 1+P1_1 \ 

«Ai— 0—Aui - / d 1 a 2 g + i Ati 2 + (p x v - tti/i) -g- + OlOil— g —- 

«2 « 2 1 « 2 1 U 2 


V 


<■-' , . < 
+a:i<J — - 5 —dit — b-2Pi——_ 

îi , 1 m, 


( n+ai 
n_ 

f Sl+/?l+l 


dil. 


) 


\ Z Z / 

En appliquant la formule de Green pour les termes en Laplacien, on obtient 

u, [fduAuf 1 - 1 

W VM!J + W(+)Àr 


,cti — 1 




a ,n 1^1 < 1_1 

-p —A Uid\l = —aiüi f V ( —; 


dT, 


66 



du 

et co m me —- = 0 (condition de Neumann), 
or] 

— Ot\CL\ 


alors 


o 


" U ai ~ 

— AuidQ 



2 «j 1 Vu, -/3 1 u" 1 - 1 m5‘ _1 Vii 2 ) ViddSi, 


donc 


" u ai ~ l f ( 

aiai—, — AuidÇ} = / —ai (ai — 1) oi—— 

'4‘ l \ «2 


n 


u ai ~ 2 u’ yi ~ 1 ' 

— |V?ii| 2 + ai/3^1 \ +l VuiVu 2 ) dfl. 

«2 ufi 


( u^^ \ 

Pour le terme fi 1 a 2 —^—^Au 2 1 <A2, on trouve 
1 Uo 1 J 




v ai \ 

f3ia 2 -J— ï Au2 dû 

V J 



■u" 1-1 u ai 

-aifi l a 2 i 1 Vui'Vu 2 + fii (fii + 1) «2 g +0 
u 2 1 U 2 l 


n ai 

1 |Vu 2 | 2 


dU. 


Alors on obtient 

u ai ~ 2 u ai ~ l 

l)ai^—|V«i| + aifii (ai + a 2 ) — 7 — r V-uiVu 2 
u 2 1 ufi 

u ai u ai 

—fii (fii + 1 ) a 2 /3 2 +2 | Vn 2 1 2 + (fi^ — cti/i) —^7 

U 2 

f ri+a!i 


-ai (ai 




\ 


V 


z/ - ai/i) -4- 
n 2 

^ai+pi-l M ri+ai n “ 1_1 

+Ml .a+ 91 _ s'+^+i + 

a 2 a 2 a 2 


dfl. 


J 


J = 



1 ? .a2-l ? .^2 f d U 3 

u f* v 2 3 4 v àt , 

O 4 


■ fi 2«3 2 «f 1 


' \ 
a 3 A « 3 - /i?i 3 + 63 -fi- + o - + a(ui - u 2 + b) ) 
, U A J 

/ rfï \ 

^ ^a 4 A « 4 - vu 4 + 64 -I 7 J 

u ? 2 . 


\ 


/ 


dfi, 


| dfi, 


/ 

n “ 2 n Q2 M a 2 +P 2 -i \ 

^+Î AU4 + (/V - «2/i) -4^ + &3«2 % 2+g2 

^ , J 4 , 

7/ “2-l 7/ Q2_1 

U/O , , x U/O 

+a 2 cr o - h a 2 a(«i — u 2 + 0 ) a 

«4 W 4 


u^ +q2 

1 ( « 2 +r 2 


dO 
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uimsciiii m iu 

i «g 2 " 1 

«2^3- g— 

^ «4 


(1Q = —«20-3 / V 


n 

<9a 3 . 

et comme —— = 0 (condition 

Oïj 

,,«2 — 1 


n 

de 


@2 

«4 


Neumann), alors 


n 


[ a2Ü3 ul 2 


Un 2 ~ f 

' 8 _ A u 3 dU = — a 2 a 3 / V 


-i ' 


n 


«f 2 


Va 3 df 2 , 


n 


Pour le terme J 

ü 


, m «2 
/ 5 2 a 4^f- î AM 4 

V < 2 


cto — 

«2 (<22 — 1) ®3 3 n 

U A 2 

cIO, on trouve 


^«2-2 ^a 2 -I 

|Vu 3 |" + « 2 ^ 20-3 g + 1 Vu 3 Vn 4 

«4 u 4 2 


n 


f . Mo 2 . 1 

, P 2 «4^tyAm 4 

( U. 4 2 


dfi : 


n 


alors J devient 

/ u° 2 ~ 2 

-«2 (CÜ 2 ~ 1) 0 3 ~4— |V« 3 | + 

M 4 2 


H «2 ^ l / a2 

■/5 2 a 2 a 4—Vu 3 Vrt 4 + P 2 (P 2 + 1) |Va 4 | 

U 4 U 4 


I 2 1 df 


J: 


n 


n“ 2_1 „ 

a 2 P 2 («3 + o 4 ) -4— T Va 3 Va 4 

«4 

rvo Cüo-I-Do— 1 


Un 2 9 Mo' 

-P 2 (.P2 + !) °4-^T2 \^ U P~ + (P 2 V ~ «2/i) -g- 

M 4 2 ^ M4 2 M4 


V 

Donc 




n a2_1 u a/2+r2 ~ 

+a 2 a( Ul -u 2 + b)^j^ - b 4 p 2 ^ +J2+1 
U 4 114 


1 

,«2 ..« 2 +P 2 - 

«3 0,3 

+ « 20 -- 


«3 , î £ 2 

C64 ce 4 


U?- 


I + J devient 


( 


ïM'» : 


O 


u ai ~ 2 a * 1 - 1 

A—— |Viii|“ + «i/3 1 (ai + a 2 ) 1 + 1 V«iV «2 

7 / 1 , 1 a!; 1 

ao-2 


_ 0 (| (c^i — 1 ) —g - 

U 2 

U a 1 2 m «2- 

-Pl {Pl + 1) « 2^-2 | Vu 2 | 2 - «2 (a 2 - 1) «3-V" 

«2 1 M4 2 

ry o — l 


,/»-2 

■Aï- lv “ 31 

,C*2 


< 2_1 

+«2/?2 (°3 + o 4 ) g + i Va 3 Va 4 - 

7/4 

, < t < 1+P1_1 , ai 

" - ' ' “ Mi 51+ft+1 + a,a 

u 0 1 a 


+ (/4 


< 2 

-Va 3 Va 4 - P 2 (/^2 + !) a 4~â^ |Va 4 |' 

u 4 

ri+ai ai 

‘zi 


U 1 U p 

V ~ " 1/0 + & 1"1 V+gi - „ 

«2 M 2 % 

,«2 «2+P2-1 

60 7 Cl/O U/O 

,z/ - a 2 /i) -g- + b 3 a 2 —„ + a 2 ci— j~ 

Ia Ua < 2 


+ {P 2 1 


uf 


-1 


V 


4*«20-('Ui 


U 4 

,,« 2-1 „«2 + /'2 

U2 + b)^-b 4 P 2 ^^ 

«4 «4 



avec 


= h +h 


h = 


n 


/ u ai ~ 2 m “ 1_1 

! —a± (ai — 1 ) a\ 1 ^ |Vwi | 2 + ai/?! (ai + 02) g i+1 

I M 2 ll 2 

v ai 

-Pi (Pi + l)a 2 ^\Vu 2 \ 2 

ti, 


U? 1- * w \ 

■Vu,iVm 2 ' 


V 


dQ, 


O 


/ 2 uT~ l 

/ -a 2 (a 2 -l)a 3 — g— |V« 3 | + a 2 ) d 2 («3 + a 4 ) g +1 Vtt 3 Va 4 

u i U T 

1 CK 2 

~p2 {P2 + 1) °4 ^+2 \ WU *\ 2 

u 4 


J 

\ 


V 


dn, 


Ji = 



U 4 

X 1 , < 1+P1_1 , 0 «? +ai «r » Jn 

(/M - ai/i) -g- + Mi— J— - b 2 Pi +/3 +1 + a i a —I diï, 

u 2 u 2 u 2 1 u 2 


! 


ai —1 


n 


J 2 


n 


( U Œ2 u °2+P2- 

(/V - a 2f i) + b 3 a 2 

U4 LL4 

u a2 ~ 1 u a2+r2 

+a 2 a(ui -u 2 + b)-^ -M 2 ^ +S2+1 

u 4 u A J 


«2+P2-1 

u 3 1 “3 

+ ot 2 a—g- 

u 4 2 

t Ct2+T2 


dil, 


où 


ù Jt = 7i + J i et j t Mt)) 


2- 


Estimation de 


Ji 


■h 



?/ «l ai+Pi-l ri+ai ?/" 1_1 \ 

(/V - ain) + Mi - b 2 Pi s ^ + g i+1 + aidû, 

u 2 u 2 u 2 u 2 J 


,“1 / > 7/ ai+Pl_1 /■ 7 / ri+Ql /■?/“ 
+ 6i«i / - Ml / + « W ~ 

J lin J lin J lin J U 

n n n 


ai —1 


-s—dQ. 

«2 
n 


fu a 

= (/V - ai fi) / 

J lin 

n 

D’après le principe de maximum il existe c 3 qui dépend de HMloo de IIMI 
que w 2 > c 8 > 0. On ainsi 


telle 


«i — l 


«i — l 
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donnant 


«i — l 

< 1_1 (uT\ “i 

dT — ° 9 1 “737 

u 2 \ u 2 J 


OU Cg = — 



01 

«1 


et en appliquant le Lemme 4.2.1 on obtient 


ai+pi —1 
U 1 

U *+« 


/ o < 1+ri i u r 
- ^ 7 /i+*i+i + c I ..01 
a 2 


0 i 


u 


alors 


&iaq 


Ql+Pl-l 

_ 

d 1+, ‘ 





où Cio = b\C. 

Donc d’après (4.2.34) et (4.2.35), on obtient 


fv ai f 7/“ 1+ri f (v ai 

J\ < (01^ - Oiifj) / + b 2 /3 1 / j3 | +Si+1 rffi + / Cio I 

«/ ’lio «/ \ Un 


+QLl<7 



I 11 J_ ' 

dü — b 2 /3i / p i+Sl+1 d£l, 

u 2 


cci+ri 


alors 





donc 


a. î 




cio ( -4- ) dD < en (cci(t)) 
«2 


01 


OÙ Cn = Cio |^| 1_ 


(4.2.34) 


(4.2.35) 


dn 


î 


dn, 
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On a 





ai 


5 


alors 




«i — l 


a i — l 

JO < C12 (aq(t)) °ù 


OÙ C12 = Cg |H|"l , 


il s’ensuit 


«i — l 

Jl < (Pi" ~ «1 II) U! (t) + C n (^1 (t)) 91 + «i(TCi2 (u>i(t)) "1 


ce qui donne 


Ji < (3/3^ - «i/i) wi(t) + ci 3 


/ ai~l \ 

I (oqfJ)) 01 + «icr(o;i(t)) a ± j 


V 


J 


Estimation de J 2 

On a 


u 


OL 2 


J 2 = (/3 2 Z/ - «2/U.) / -4-dll + 6 3 «2 




«2+P2-1 






/32+<?2 


-dfl + «2<T 




02 — 1 




u A 


U 


02 — 1 


U 


02+r 2 


+« 2 a / («1 - u 2 + 6)—dfi - 6 4 /3 2 / g JH, 
d tp, 2 J U A 


donc 




02 


J2 ^ (/3 2 Z2 — a 2lP) I ~~^dn + &3«2 


U 


«2+P2-1 


U 


■u 


P2+12 


-dU + « 2 (j 


U 


02 — 1 


■u 


&2 


-dfl + 


02-1 


« 2 a / («i + b)^-dn - 6 4 /3 2 / ^ +S2+ i ^, 

1/4 


t/ ^ «/ 

n n 


02+r 2 
3_ 

2 " 


et comme nous avons montré dans le chapitre précédent que les solutions du système 
composé de deux premières équations du système (4.2.27), existent globalement et 
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bornées (3 H > 0, telle que a\. «2 < H), donc on peut écrire 


J 2 < {/3 2 v - «2/-0 / -4-dD + 6 3 a 2 


/ OL2 /' ...«2+P2-1 é?/ 02-1 

^ + &3«2 / -^^dîî + « 2 a / + 

Ma J U A 2 J U A 2 

n n n 


a 2 a(iJ + 6) / -^j^dü - b A ( 3 2 / ^ +82+1 dfl. 


D’après le principe du maximum il existe c 44 qui dépend de ||IIoo ,? Il^2II00 > Il ^3 II 00 
de ||</? 4 || telle que n 4 > C14 > 0, donc 


a 2 - 1 


Mg 2 \ «2 


- r 2 < 


«2 - 1 


uf \ « 2 


" Wf 2 


« 2 -l 

u" 2 " 1 /< 2 \ « 2 

if sci5 U 


OU C15 — I — 
\Cl 4 


D’après le Lemme 4.2.1 


02+P2 — 1 02+V2 

l _A _ < R _ 

n .02+<12 — „,P2+S2 + 1 


Q2+P2-1 


&3Q;2 é 2 +® - 6 3/?2 ^ +S2 + i + & 3 C 


- hP 2 £+*2+1 + Cl 6 


OÙ Ci 6 = 63C. 


Donc 


r u “ 2 r m “2 +r- 2 /• (u a2 \ 

J 2 < (/V - « 2 /x) / -\dü + b A (d 2 / / 3 2 3 +Æ2+1 dQ + / ci 6 ^ dD 

J U A J U A J \ U A J 

n n n x 7 

«2-1 

r ( ,, a2 \ a, r v a 2+r2 

+«2 (cr + a(H + 6)) / Ci 5 I J dfl — b A f3 2 / / g 2+S2+1 dD, 

u V“ 4 / n “ 4 
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conduisons à 


,n , f. ( «r 


02 


J-2 < (/^ - «2 Ai) / -Jj-dü + / Ci6 - 5 - dH 






a 2 -l 


+«2 (o - + a(H + 6)) I C15 [ —g- 

«4 


W? 2 \ « 2 


s' ‘ da 


L’inégalité de Hôlder entraine 

#2 


MS *d/ S 1 " 

n ^ 4 / Vu ' / \n 


02 / \ 1 — 02 
(cie) 1 ^ dfi 


donc 


/ n \ 

( u 2 \ 

Cl6 -f- dn < c 17 (u 2 (t)) 02 , OÙ Ci7 = Ci 6 |fi| 


1 — 02 


U 


A 2 J 


On a 


« 2 -l 


HS) “ 2 


a 2 -l 
«2 


_j_ 

CK2 


(ci 5 r dn 


alors 


« 2 -l 


< 2 \ «2 


« 2 -l 


Ci5 I -g- 1 dO < Ci 8 (u; 2 (t)) "2 , 

«4 2 


OÙ Ci 8 = C15 |0|“2 . 


On obtient 


«2-1 


d 2 < {P 2 v - «2/i) U 2 {t) + Ci7 (u 2 (t)) 92 + «2 {cr + a(H + b)) ci 8 (cu 2 (t)) " 2 , 


ce qui donne 


d 2 < (^ 2 ^ - « 2 ^ 2 (t) + Cig 


/ 


(u;2(d)) 02 + «2 (<7 + a(H + 6)) (u; 2 (t)) a2 


V 


«2 - 1 \ 

/ 
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Revenons aux intégrales I\ et I 2 

ai—2 


( U ai ~ Z M ai_1 \ 

/ R, (01 - 1) », l Q \Vui\ 2 + a 1 /3 1 (ai +a 2 ) l +1 'Vu 1 'Vu 2 \ 

Uo 1 Mi 1 


h = 


U. 

-Pl (Pl + 1) a 2 p i+2 |V«2 


v \y\ 1 -v 

\ U. 

( — OL\ («i — 1) CI1U2 |V«i| 2 

+a 1 p i (ai + a 2 ) UiU 2 Vwi Vm 2 

y -/?i (/?i +1) a 2 «? |Vm 2 | j y 

f Ui 1 ' 2 

~ I Q l /q+2 6 ^’ 


CÙC 


■u" 1 2 

/^i+2 

«2 


dx 


u. 


n 


ou 


/ 


Qi — 


«i («i — 1) ai«| |Vui| 


\ 


—aiP 1 (01 + o 2 ) UiU 2 WU\Wu 2 

\ +Pi(Pi + l)a 2 u\\Vu 2 \ 2 J 


Qi est une forme quadratique par rapport à: u 2 V//] , U\ V// 2 qu’on peut écrire sous 
la forme matricielle suivante: 


Q 1 — 


^ «i (ai — 1) ai —-ai/3i (ai + a 2 )^ 
1 2 




\--a 1 p i (ai + a 2 ) /5 1 (/3 1 + l)a 2 y 

Q, est définie positive si, et seulement si, tous les déterminants principaux successifs 
Ai, A 2 de sa matrice des coefficients, sont positifs, 
où 

Ai = oq (oq — 1) cp, 

Ai > 0 d’après (4.2.31). 


A 2 


ol 1 (et 1 — 1) a x — oti/3 1 (ai -f- a 2 ) 

1 2 
— 2 a lPl ( a l + ÛS 2 ) /?! (/?! + 1) Û2 


a\f3\a\a 2 


(«î — 1) (/?i + 1) / ai + a 2 \ 


«i 




v2\/ a i°2 / 
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Nous avons A 2 > 0 d’après (4.2.32). 
Donc 




ou 


/ / Mo 

1 -«2 («2 - 1 ) a 3 


« 2-2 


h = 


U 




| Vn 3 | 2 + a 2 /3 2 (a 3 + o 4 ) 


u 


02—1 


/?2 + l 


Vu 3 V?i 4 


\ 


U 


~@2 (/^2 + 1 ) a 4 


.«2 


^ 2+2 


|V«4 


/ - 


«2 («2-1) a 3 «4 |V« 3 | 




+« 2/^2 (°3 + ) U 3 U4VU 3 VU4 

^ ~@2 (^2 + !) a 4«i |V?i 4 | 2 y 


U 


02 — 2 


/? 2+2 


<ir 


u 


U 


Ol 2 — 2 


^2 02+2 °^5 


U 


Q 2 — 


^ «2 («2 - 1 ) 0 3 « 4 | Vîi 3 | 2 ^ 

—a 2 f3 2 ( a 3 + a 4) M 3 M 4 Vm 3 Vm 4 

y +/3 2 (/3 2 + l)a 4 ^ |Vn 4 | 2 y 


dfi, 


Q 2 est- aussi une forme quadratique par rapport à: «4 |Vu 3 | , u 3 | V 1 x 4 1 ,on a 

Q 2 = 


^ «2 (&2 ~ 1 ) a 3 ~^ a 2/^2(03 + 04)^ 

1 


-a 2 /3 2 ( a 3 + a 4 ) (3 2 (P 2 + 1) a 4 J 




Q -2 est définie positive si, et seulement si, tous les déterminants principaux successifs 
A 3 , A 4 de sa matrice des coefficients, sont positifs, 
on a 

A 3 = «2 (Cl 2 — 1) 0 3 , 
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A 3 > 0 d’après (4.2.31). 


A 4 


«2 {pt 2 — 1) a 3 — ~Ol-2^2 ( a 3 + 0 4 ) 

1 

— 2 a 2^2 (°3 + 04 ) P2 (P2 + 1) a 4 


^2/^ 2 ^3 ^'4 


( 1^2 — 1) (/3 2 + 1) 
«2 /?2 


/ a 3 + a 4 \ 
\2^y o 3 o 4 / 


la condition (4.2.32) nous donne la positivité de A 4 . 
Donc 


— (u 2 (t)) < J2 ■ 


? 


Enfin on a 


| (wW) = | (wiW) + | M()) 
< J\ + J2 

/ 


< (/? 4 z/ - ai/r)o;i(t) + c 43 

/ 


(wi(t)) 01 + «îcr (cni(t)) “i 


V 


(o^)) 02 + «2 (cr + a(H + 6)) (cu 2 (t)) a2 


«i - 1\ 

/ 


V 


+ (P 2 v - a 2 n)u 2 (t) + C19 
d’après (4.2.31) on aura 

w'iW < cis^i 1 ^) + ci 3 ai<r(a;i (t))^ü“ 


a 2 - 1\ 

/ 


(4.2.36) 


Par ailleurs, on obtient par (4.2.31) l’inégalité 

oj' 2 {t) < c 19 u e 2 2 (t) + Ci 9 a 2 (cr + a(H + b )) (o; 2 (£)) “ 2 • (4.2.37) 

Par une intégration simple de (4.2.36) et de (4.2.37) dans [0, T*] en utilisant le Lemme 
4.1.4 on se rend compte que u(t) est uniformément bornée; i.e. uj(t) < q 4 ,où q 4 dépend 

de H^ilL , Moo » Hélico et de ML- ■ 

Preuve du Corollaire 4.2.1. En appliquant le Théorème 4.2.1 sur les termes 
de réaction, qui sont dans L°° ((0, T max ), L m (O)), pour tout m > A(A±i) ^ moyennant 
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les résultats de Haraux-Kirane [16] et de la Proposition 48.4 de [43], on conclut que la 
solution du système (4.2.27) est globale et uniformément bornée dans fl x (0, oo). Voir 
aussi la preuve du Théorème 3.3.1 pour plus de détails. ■ 

Remarque 4.2.1 Si on pose: 


CL\ 

d Ul : 

— d U2 , O3 

du 3 ) 

a 4 = d U41 a = p 0 p, a = ky 3 , 

h 

— c 0 p, 

b 2 = dp’i 

h = c 0 p, 

b 4 = dp', b = floor, 

Pi 

= 2, qi 

= 1, n = 2, 

si = 0, 

P2 = 2, q 2 = 1, r 2 = 2, s 2 = 0, 


on obtient le système (4-2.26), où ses constantes réalisent les conditions (4-2.30) de 
l’existence globale des solutions en temps. 
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Chapitre 5 


Etude de l’explosion en temps fini 
des solutions d’un Système Général 
de Gierer-Meinhardt à m 
composantes 


Comme nous avons vu dans le Chapitre 4 les solutions de système de Gierer-Meinhardt 
avec trois composantes explosent en temps fini sous des conditions sur les données 
initiales de système EDO associé, alors nous allons généraliser l’étude de l’explosion 
en temps fini pour le système EDO associé au système générale de Réaction-Diffusion 
de type Gierer-Meinhardt de m équation ou en terme biologique; un modèle composé 
d’un Activateur et de (m — 1) Inhibiteurs avec des concentrations (uj), j = 1, 

Soit le système 


dui u^ 11 

—- = CtiAtii - Mi + -s; - h <y 

ot 

drij 

— = ajAuj - b jUj + m 

dt nA- 

i =2 * 


m 

n vA li 

i =2 * 

uT 


,x £ Q ,t > 0 


J = 2, ...,m 


avec les conditions aux bords de Neumann 


(J LL ' 

—- — 0, xG(9f2,t>0, pour j — 1 ,m 

ori 


(5.0.1) 


(5.0.2) 
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et les conditions initiales 


Uj (. x , 0) = (pj(x) >0, x G fl, pour tout j = 1, m. (5.0.3) 


Comme nous avons énoncé précédemment, fl est un domaine borné de classe C 1 dans 
M. N , de frontière T = d fl. 

Supposons que a > 0. aj,bj,Pji des indices non négatifs pour tous: i,j = 1 
avec 


Pu - 1 ^ . 

- - -- > mm 

max [pj i) 

j=2,m 


min (pji) 

j,i=2,m 


.np« 


m ’ 

max (pu) + Il pu 

i=2,m 


(5.0.4) 


L’explosion des solutions se réalise s’il existe un temps T ma x < oo, tel que les 
solutions sont bien définies pour tout 0 < t < T max , où 


m 


5Zlk(i,.)ll 0O = °o- 

T'/ J- max 

î= 1 


Lemme 5.0.2 Soit (ni (t, .),..., u m (t,.)) la solution de (5.0.1)-(5.0.3). Alors pour tout 
(■ t , x ) dans (0, T max ) x fl, on a 


u l (t, x) > e bit min (</q(x)) > 0 , i = 1 ,..., m 


(5.0.5) 


Théorème 5.0.2 Supposons la condition (5.0.4), alors il existe un espace-indépendant, 
des conditions initiales i = 1 ,..., m, tel que la solution (ni,..., u m ) = (ui(t ),..., u m (t)) du 
problème (5.0.1)-(5.0.3) explose en temps fini T max . 


Preuve. Considérons Pespace-indépendant des solutions de (5.0.1)-(5.0.3), i.e. les 
solutions de système EDO associé sans diffusion. Pour des conditions homogènes spa¬ 
tiales > 1, i — 1 supposons que T max (p 11 ..., <p m ) > 1 (pour obtenir une 

contradiction). Dans ce qui suit, toutes les constantes positives C) sont indépendantes 
de <p t , i = 1,..., m. 

™ n Q1 

Pour a* > 0 , i — 1,..., m, soit À = aq/q — \ . «A et L ( t ) = —. 

i =2 n ufi 

i=2 1 

On a 


ü(t) = ai ^— 

.w* 



m 

E 

i=2 


U 


a,- 


a 1 

1 


U. 


<u+1 


m 

n u 

3 =2 


i 



5 
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se iï 


Oi,\ 


utilisant les équations de (5.0.1) sans diffusion, on obtient 


/ 


L\t) = 


a i rn 

7 ^1 \—> , 

— ' blOil— -h 2__j bi a i 


m 

n u ai 

i =2 


i= 2 


pn+ai-1 
(Xi V > 

+ a lxs- 


U y 1 
m 

n u? 1 

i =2 


U 


Pjl+Ol 


n « 

i =2 


,Pli+oti 


aj- 


j =2 U 


+ cra i“^— 

n 

i=2 * 


aj+Pjj+i A Pji+a, 


n vt 

i =2 

i¥=j 


donc 


?/* 

L'(f) + AL(*) = 


.Pii+ai-l 


E 


a i 


Pi i+ai 

U 1 


n w ?n+«i ’ Pli II 

i=2 * J 


n < 

î=2 * 


u? 1 " 1 

+ CTQ! iisr 


m 

n<* 

i=2 * 


Nous envisageons deux cas séparés. 

Cas 1: 


(5.0.6) 


pu - 1 > max (pji). 

j=2,m 

Appliquant (5.0.6) avec oq = 1. prenant cq suffisamment grand, i = 2, m , et utilisant 
(5.0.5), (pj(x) > 1, j — 2, m, et l’inégalité de Young on aura pour tout t G [0,1] 


«f 1 


/ 


u 


,Pli+Oli 


OL\ 

m 


U 


t pj l+l 


\ 


Pll/Pjl+1 


va«r'""y 


Pjj + l+Qi 


m , 

n U; 13 > aq 

Î=2 J 


u 


Pli+ai 

1 


U 


777 / 3 

a j +Pjj +1 JJ yPji+ a i 
i =2 * 


Ci , j = 2, m 


ou 


hij — (Pji Y 1 + °o) 


Pli 


Pjl + 1 


çij — «j > 0 , j = 2, m, i = 2, m 


Cj- = j 

V m / J 

D’autre part en utilisant (5.0.5) on aura 


= 2 , m. 


«i -u^ 11 

m nuf i+ai 

î=2 * 


«i 

m 


/ \ 


Pli 


^1 


777 , 

v,a<v 


n «;■ > C m+ , 

1=2 


( V 11 

Ui 
m 
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ou 


di = (pu ~ 1) OLi - pu > 0 „ i = 2, m. 


Alors on obtient 


aq- 


n u 

i =2 


pu ™ 

n- • -L/-v • ‘ — 


U 


Pji+ai 

1 


Pli+Ctj 


i=2 U 


a j+Pjj + 1 TT n ,Pji+ a i 


+ a 


ra+1 


n< 

i=2 * 

*A? 


/ V“ 

Ui 
m 

V,5<7 


-E c / 

i=2 


de (5.0.6), nous avons 


/ \ 


Pli 


L'if) + AL(t) > C m _|_i 


«i 


v.a<7 


E°j 

3 =2 


ou encore 


où C = XI Cj. 

3 =2 


P(t) + AL(t) > C m+ i (L(f)) pn - C 


Pour L(0) suffisamment grand, u i explose avant t = 1 ; ce qui contredit la suppo¬ 
sition T max (<p x ,..., v? m ) > 1. 

Cas 2: 


Pu - 1 < max (pj i) , 

j=2,m 


min_ (pji) - max (pi*) + lIp H 

j,i=2,m i=2,m ^—2 


(Pu - !) 

Affirmons qu’il existe des constantes C m+ 2 , C m+ 3 > 0 telles que, si 


> Ilpijmax (pji) ■ 

*—2 j=2,m 


max_(Pj 1 ) —pi 1 +1 m min_(pij)-gij ) 

/ > r ïïi/i p =2,m ' 

Pl — '- / m+2 Pi 

3 = 2 


(5.0.7) 


alors 


max (p,i)-pn+l ) min (p,,) // 1; ) 

\j=2,m I s-sj t-t \i=2,m I _ , . . 

u i■ 7 > C m+3 lia) 7 , pour 0 < t < 1. 

7=2 ù 


(5.0.8) 


Pour prouver cette affirmation, prenons 0 (t) = e Af L (t) et appliquons (5.0.6) avec 
«i = max (pji) — pu + 1 > 0, 

j=2,m 
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oij = min (Pij) - qij , j = 2 ,m. 

2=2,m 

Alors, pour tout t G [0,1], 


donc 


et ainsi 


$ (t) < 0 ==► U ( t ) + AL (t) < 0, 


max (pn)+ai 

7=2,m 


U 7 


min (p ? ;;)+g; ) Il U 


m | 3 

\ z=2, m 


n u 

3 =2 J 


> 


M' 


pn-l+ai 

1 




j=2 3 


vi=2 


a l I / , 'G? I 3 


*(f)>e-W^>rW«,A 


-1 


Ct; 


=:C, 


m+2- 


n v 

j=2 3 


0=2 


Par conséquent, on a 0 (t) > (7 m+ 2 dans [0,1], et l’affirmation est prouvée en prenant 

Cm+3 — e ' A ^m+2- 


Maintenant, supposons (5.0.7). Utilisant la première équation dans (5.0.1) et 
(5.0.8), on déduit 




^ 1^1 ^ m ^ ^ 772 + 4^1 


/ n >»Ë_(pj 1 )-p 11 +i 

Î2 1=2 .m' ' 

î* 11 ” P Plj 7 7 Y 7 V 

J — 2 mm Pii “ max Pi ; 

\ j,i=2,m y ' j=2,m v ' 


n« 

i=2 3 


(pu) 


= C' m+ 4'«1, 0 < £ < 1 


ou 


2,P 

( - / m+4 l -'m+3 


(pM- max_( Plj ) 
x 7 j=2,m v 7 


0 = pu - np y 


= 1 + 


max (pji) - pu + 1 

j=2,m 

j= 2 ^ 1J min (pjj) — max (pij) 

j,i=2,m j=2,m 


(fin - 1) 

m 

min _(pji) - max (pi;) + IIpi, 

m 

~ n pij max (p fi ) 


N 

II 

<-> 

II 

.11 

3~ 2 j=2,m 


min (pji) — nrax (pij) 

j,i=2,m j=2,m 


> 1. 


En prenant cp, suffisamment grand, on obtient l’explosion de » i avant t = 1; ce qui 
contredit la supposition T max >1. ■ 
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Conclusion 


Dans ce travail nous avons établi l’existence globale des solutions d’un certain sys¬ 
tème de Réaction-Diffusion de type Gierer-Meinhardt composé de trois équations mod¬ 
élise un phénomène biologique appelé Phyllotaxie; 


du . u pi 

— = ai Au — a i u H- 7 -— 

dt v qi ( w ri + c) 


. dv . u p2 

' ttt = a 2 Av - p 2 v H-. 

dt v&w r 2 


dw u P3 

= a 3 Aiu — u,nw H-, 

dt ^ v q zw r 3 ’ 


+ a i 


sous la condition 


( q i rq \ ( r i q\ 

0 < pi - 1 < max < p 2 mm -—-, —, 1 ,p 3 mm -—-, —, 1 

l \Q 2 + 1 r 2 J \r 3 11 g 3 

D’autre part, nous avons démontré l’explosion en temps fini des solutions de son 
système EDO associé à condition 


Pi 


max (p 2 ,p:i) 


> mm 




,1 


min (ç 2 , 93 , r 2 , r 3 ) - max (<?i, ry) + qpri : 

Ce dernier résultat d’explosion peut être généralisé pour le système de m équations 
qui est un modèle composé d’un activateur et de (m — 1) inhibiteurs; 

( dui 


= GqA U\ — biUi + 


u , 


dt 

dUn . . U,1 

— = ajAuj - b jUj + — 

at n u 

i=2 


+ 


Pli 


n u: 

i =2 " 
Pi 1 

up 


Pji 


J = 2,..., m 


avec 


P il 


max (pji) 

j=2,m 


> mm 


m 


min (p^) - max (p u ) + n p u 

j,i=2,m i=2,m ^—2 
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